Astronomy & Astrophysics manuscript no. 1643 


© ESO 2009 


March 11,2009 





0^ 
O 
O 
(N 



O 



> 

o 

rn 
O 

a^ 
o 



>< 



Sub-surface convection zones in hot massive stars and their 

observable consequences 

M. Cantiello^ N. Langer^-^, I. Brott\ A. de Koter^'^, S. N. Shore"*, J. S. Vink^, A. Voegler*, D. J. Lennon^ and 

S.-C. Yoon^ 

' Astronomical Institute, Utrecht University, Princetonplein 5, 3584 CC, Utrecht, The Netherlands 
Argelander-lnstitut fiir Astronomie der Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany 
Astronomical Institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 SJ, Amsterdam, The Netherlands 
Dipartmento di Fisica "Enrico Fermi", Universita di Pisa, via Buonarroti 2, Pisa 56127 and INFN - Sezione di Pisa, Italy 
Armagh Observatory, College Hill, Armagh, BT61 9DG, Northern Ireland (UK) 
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

Department of Astronomy & Astrophysics, University of California, Santa Cruz, High Street, Santa Cruz, CA 95064, USA 



Received 12 January 2009 / Accepted 03 March 2009 

ABSTRACT 

Context. We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with 
iron and helium ionization. 

Aims. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront 
our results with observations of OB stars. 

Methods. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the 
mixing length theory is used to characterize the envelope convection zones. 

Results. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial 
metallicity. It is absent for luminosities below about 10'^ Lg, 10^ ' Lq, and 10*-^ Lq for the Galaxy, LMC and SMC, respectively. We 
map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of 
observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. 
Conclusions. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities 
argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere 
of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same 
mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption 
components in ultraviolet absorption lines. 

Key words. Convection, Stars: atmospheres. Stars: magnetic fields. Stars: oscillations. Stars: mass-loss 

1 . Introduction While the envelope convection zones may, at first glance, be 

negligible for the internal evolution of hot massive stars, they 
Massive stars, in a general sense, have convectiv e cores and ra- may cause observable phenomena at the stellar surface. The rea- 
^ diative envelopes dKippenhahn & Weiger3ll990h . pie introduc- son is that the zones are located very close to the photosphere 
tion o f the so called "iron peak" in stellar opacities jlglesias et alj for some mass interval (see below). Here, we will discuss which 
" " ' 119921) led, however, to the prediction of a small convection zone observed features in hot stars might be produced by these near 
in the envelope of sufficiently luminous massive main sequence surface convection zones. In particular, we examine whether 
models (Stothers & Chin 1993). It is often accompanied by an a link exists between these convective regions and observable 
even smaller convection zone which originates from an opacity small scale velocity fields at the stellar surface and in the stel- 
peak associated with partial helium ionization. These two con- lar wind, "microturbulence". A similar ide a has been used to 
vection zones comprise almost negligible amount of mass. The explain microturbulence in low mass stars jEdmundsll 19781) . in 
realit y of the iron opacity bump, as predicted by various groups which deeper envelope convection zones reach the photosphere. 
(e.g.. |lglesias et al. 1992; Badnell et al. 2005), is unambiguous. While Edmunds ( 1978) concludes that the same mechanism can- 
It is most obvious in the field of stellar pulsations. Only the not explain microturbulent velocities in O and B stars, the iron- 
inclusion of this feature allows an agreement of observed and peak induced sub-photospheric convection zones in these stars 
predicted instabilit y regimes in the HR diagrani, fro m the white had not yet been discovered. We demonstrate in this paper that 
dwarf regime (e.g. ISaiolll993t ICharpinet e t al.' 1997), for main these convection zones may not only cause motions which are 
sequence stars (e.g., fi Cephei stars; see Deng & Xiong 200 L observable, but possibly even directly affect the evolution: First, 



and r eferences therein), and up to hot supergiants dSaio et alj we discuss how photospheric velocity fields may affect the struc- 



120061) . ture of massive star winds by inducing clumping at the base of 

the wind and thereby affecting the stellar mass-loss. And second, 
we argue that the near surface convection zones may generate 
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affect the stellar wind mass-loss and, more significantly, the as- 
sociated stellar angular momentum loss. 

We construct grids of massive main sequence star models, 
for various metallicities, that allow us to predict the occurrence 
and properties of sub-surface convection zones as function of the 
stellar parameters (Sect. |3]l. We then compare the model predic- 
tions with observed stellar properties, e.g., empirically derived 
microturbulent velocities and observations of wind clumping in 
hot massive stars (Sect.|4|i. 

2. Method 

Our stellar models are calculated with a hydrodynamic stellar 
evolution code. This code can calculate the effect of rotation on 
the stellar structure, rotationally induced chemical mixing, and 
the transport of angular momentum by magnetic torques (see 
iPetrovic et al.ll20()5t lYoon et al.ll2006i and references therein). 
Compositional mixing is treated as a diffusive process. The rate 
of change of a nuclear species of mass fraction X,- is calculated 
as 



\ dt I 



[dm 



dt 



(1) 



where D is the diffusion coefficient constructed from the sum 
of individu al diffusion coeffi cients for the range of mixing pro- 
cesses (see iHeger et al.l 12000 , and references therein). The sec- 
ond term on the right hand side is the schematic symbol to stand 
for all nuclear reactions. The contributions to the diffusion co- 
efficient are convection, semiconvection and thermohaline mix- 
ing. For rotating models also the contributions from rotation- 
ally induced mixing and magnetic diffusion are computed. The 
trans port of angular momentum is also treated as a diffusive pro- 
cess (lEndal & Sofialll978t iPinsonneault et al.ll 19891; iHeger et all 
120001) . 

The Ledoux criterion is used to determine which regions of 
the star are unstable to convection: 



ad 



• V + ^ V„ 



< 



(2) 



(e.g.. iKippenhahn & Weig"er3ll990h where Vad is the adiabatic 
temperature gradient and V^, is the gradient in the mean molecu- 
lar weight. The diffusion coefficient, D, in convective regions is 
approximated with 



D = —aHpVc 



(3) 



where Hp is the pressure scale height, Uc is the convective veloc- 
ity, and a the mixing length parameter. We fix a - 1.5, which 
results from evolutionary tracks of the Sun (e.g. ' Abbett et al.l 
119971; iLudwig et al1ll999l) : a sensitivity study of the a depen- 
dence of our scenario will be presented in future work. The con- 
vective v elocity, Vc, is calculated using the mixing length theory 
dBohm-Vi tense 1958) (MLT hereafter) and the convective con- 
tribution to the diffusion coefficient becomes: 
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— g/3(l -/3)Vad(V,ad-Vad) 



1/3 



(4) 



where x is the opacity, p is the density, /3 is the ratio of gas pres- 
sure to total pressure, g is the local gravitational acceleration, 
and c is the speed of light. Here, Vrad and Vad are the radiative 
and adiabatic gradients, respectively. 

We use the solar composition proposed by lAsplund et al.l 
(l2005h . The opacities in our code are extracted from the OPAL 
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Fig. 1. Opacity in the interior of 60 M© zero age main sequence 
stars of various metallicities (see legend) as a function of temper- 
ature, from the surface up to a temperature of 10^ K. The differ- 
ent colors refer to different metallicities, as shown in the legend. 



tables ("iglesia s & Rogersll99 6'). Fig.[T]shows the opacity coeffi- 
cient as function of temperature in our 60 Mq models for various 
metallicities. The peaks at log T ^ 4.7 and log T ^ 5.3 are 
caused by helium and iron, respectively. The peak at log T ^ 
6.2 - 6.3 is caused by carbon, oxygen and iron. 

We use the m etallicity dependent mass-loss predictions of 
IVinket al.l(l2001h . 

2. 1 . The helium convection zone 

In the very weak helium convection zone, radiative diffusion 
is the dominant energy transport mechanism, which may have 
consequences for the development of convection. In fact, in vis- 
cous fluids the Ledoux-criterion is not strictly correct, since it 
ignores any dissipative effect on the evolution of a perturbation. 
A more accurate criterion can be expressed in terms of the non- 
dimensional Rayleigh number, Ra which for compressible, strat- 
ified convection, is 



Ra:^ 



(V - V,,)L'g 



(5) 



Here L is the thickness of the convective layer, and k and v are, 
respectively, the thermal diffusivity and the kinematic (molecu- 
lar) viscosity (e.g. lShorelll992l p. 328). 

For convection to develop, Ra must exceed some critical 
value, Ra^. The estimate of Ra in the helium convective region 
depends on the choice o f the viscosity coefficient. For the Spitzer 
formula ( ISpitzeiiri962h .Ra > Rac, and the region can be con- 
sidered convective. In contra st, for the radiative viscosity (e.g, 
'Ki ppenhahn & Weigertll 19901 p. 445), Ra < Ra^. There is an ad- 
ditional uncertainty in these estimates since the expressions for 
the radiative transport coefficients in our models are strictly cor- 
rect only in the diffusion limit. Likewise, the value of the heat 
capacity Cp can vary by an order of magnitude depending on 
whether the radiative energy reservoir aT"^ is coupled to the in- 
ternal energy of the gas or not. Since the helium convection zone 
occurs very close to the surface in our models, these additional 
uncertainties could be relevant. 

Ideally, the properties of the helium convection zone could 
be studied through multi-dimensional hydrodynamic calcula- 
tions. However, the large thermal diffusivity poses a formidable 
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computational challenge since it makes the problem numerically 
stiff: the diffusive timescale is much shorter than the dynamical 
one, which leads to very short time steps if an explicit solver is 
used (unfortunately, most codes used for compressible convec- 
tion are explicit). Any simulation would have only limited value 
unless it includes a sufficiently realistic treatment of the coupling 
between plasma and radiation. 

In the presence of strong wind mass-loss, another consid- 
eration related to the He convective zone becomes important, 
due to the fact that it comprises only a tiny amount of mass. 
Convection can set in only if the turnover time Ttum - Hp/fc 
is shorter than the time scale for which convection is predicted 
to prevail at a fixed Lagrangian mass shell inside the convec- 
tion zone, Tconv, which is Tconv - AMconv/^- We find a critical 
mass-loss rate M ~ 10"^ Moyr"', above which convection has 
no time to develop in the helium region, since the wind is re- 
moving an amount of mass equivalent to the mass of the con- 
vection zone before a convective eddy can turn over (see Tab.[T]). 
For a metallicity Z=0.02, stars above 40 Mq cannot develop the 
He convection zone, and in a 20 M© such a layer is convec- 
tive only for 10 - 100 turnovers before convection moves to a 
lower mass coordinate. None of these concerns is significant for 
the iron convection zone (FeCZ hereafter), where convection is 
always fully developed. Moreover the convective velocities for 
the FeCZ are always found to be much higher than those in the 
helium convection zones. We disregard the occurrence of the he- 
lium convection zones unless it is explicitly mentioned. 



3. Results 

We calculated a grid of non-rotating stellar evolution sequences 
for initial masses between 5Mq and 100 Mq, at metallicities 
of Z=0.02, Z=0.008 and Z=0.004, roughly corresponding to 
the Galaxy, the LMC and the SMC, respectively. Additionally, 
we computed several models at lower metallicity. Since rapid 
rotation can change the properties of sub-surface convection 
dMaeder et alj l2008l) . we calculated a few rotating models to 
evaluate the effects of rotation on our results. These effects are 
discussed in Section lTTI 

Figures |2] and [3] show the evolution of the radial extent and 
location of the sub-surface convection zones in 20 M© and 60 Mq 
models during the main sequence phase. 

As outlined above, the He opacity bump at around 
log T ^ 4.7 is responsible for a convective zone which occurs 
close to the stellar surface and is very inefficient: only a very 
small fraction of the heat flux is transported by bulk motions in 
this region. The upper boundary is typically found at an optical 
depth in the range 2 < r < 10, where t is the Rosseland mean 
optical depth. Below this convective zone, the Fe opacity bump 
at around log T ^ 5.3 is associated with a more efficient and 
extended convective region. 

The radial extent of the FeCZ is quite substantial, i.e. a sig- 
nificant fraction of one solar radius, which corresponds typically 
to 2 - 10 pressure scale heights, comprising a mass on the order 
of 10"* Mq to 10"^ Mo, while the amount of mass between the 
top of the FeCZ and the stellar surface is around several times 
10"^ Mq (cf. Table 1). In the 20 Mq model the upper border of 
the FeCZ is located at t ^ 140 on the ZAMS, and at r a; 370 
on the cool side of the main sequence band. In the 60 M© model 
the upper border at ZAMS is located at t » 15, reaching t x 260 
during the late main sequence evolution. Convective velocities 
predicted by the MLT are on the order of 10s of kms where 
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Fig. 2. Evolution of the radial extent of the helium and iron 
convective regions (hatched) as function of time, from the zero 
age main sequence to roughly the end of core hydrogen burn- 
ing, for a 20 Mo star The top of the plot represents the stel- 
lar surface. Only the upper 1 Rq of the star is shown in the 
plot, while the stellar radius itself increases during the evolution. 
Upper panel: The star has a metallicity of Z=0.02, and its effec- 
tive temperature decreases from 35 000 K to 25 000 K during the 
main sequence phase. Lower panel: The star has a metallicity of 
Z=0.004, and its effective temperature decreases from 37 000 K 
to 27 000 K during the main sequence phase. The extent of the 
convection zones is smaller than in the case shown above, and 
the iron zone is absent for the first 2.5 million years. 

more extended zones achieve higher velocities. For a quantita- 
tive analysis, we define an average convective velocity 

1 

<^'>-=-zr ''-W^^'" (6) 

where Rc is the upper boundary of the convective zone, and 
where we set a - 1.5. 

From Figures|2]and[3] three trends for the extent of the sub- 
surface convection zones are noticeable. First, with increasing 
time during the main sequence evolution, these zones become 
more extended, and are located deeper inside the stellar enve- 
lope. This is because the stellar envelope expands, and becomes 
cooler, while the temperature of the opacity peak remains nearly 
constant. In our 20 Mo model at Z=0.02, the mass of the He con- 
vective zone increases from about 10"^ Mo to 2 x 10"^ Mq, and 
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Table 1. Properties of the envelope convection zones in our 20 and 60 M© models of solar metallicity. These are the same models 
shown in the top panel of Fig.|2]and in Fig. [3] The values in the table refer to t=6.41xl0^ for the 20 Mq model and t=2. 37x10'' for 
the 60 Mo model. 
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Mo 


Mo 




days 


days 


Moyr-' 
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He 


0.025 


0.08 


7.6 X 10^'' 


1.9 X lO^" 


1.8 X 10' 


2.5 


38 


7.3 X 10-"* 


20 


Fe 


0.08 


2.40 


3.6 X 10-'^ 


5.8 X 10-' 


1.8 X 10" 


0.25 


18250 


7.3 X 10-** 


60 


Fe 


0.24 


2.25 


1.6 X 10^^ 


9.8 X 10 ' 


8.5 X 10' 


0.83 


1570 


3.7 X 10-* 



" Mass contained in tlie convective region. 

* Mass in the radiative layer between the stellar surface and the upper boundary of the convective zone. 
^ Expected number of convective cells, A'ceiis '■= 

Convective turnover time, T,„n, := Hp/{v^). 
" Time that a piece of stellar material spends inside a convective region, T^onv := lS.M^aml^- 




1x10^ 2x10^ 3x10^ 
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Fig. 3. Same as Fig.|2] but for a 60 Mq star at Z=0.02. Note the 
different vertical scale, spanning the upper 4 Rq of the star. The 
effective temperature decreases from 48 000 K to 1 8 000 K dur- 
ing the main sequence phase. 

that of the FeCZ is growing from 2 x 10"'' Mq to 10-"* Mq. For 
sufficiently hot models, the helium convection zones can even 
vanish (Fig.|2] lower panel). Second, comparing the 20 M© and 
the 60 Mq model at Z=0.02 demonstrates that the FeCZ becomes 
more prominent for higher luminosity. This is because the opac- 
ity is not substantially changing among main sequence models 
at the same metallicity, such that a higher luminosity renders a 
larger portion of the envelope convectively unstable (both in ra- 
dius and mass fraction). Our models show that the FeCZ dis- 
appears below a threshold luminosity of about lO'^L© on the 
ZAMS at solar metallicity. Third, comparing the two 20 M© 
models in Fig.|2]shows that the extent of the FeCZ, and its pres- 
ence, depends on the metallicity. We find that for Z=0.001, it 
is completely absent below 40 M©, and at Z=0. 00001 it does 
not occur for M < 60 M©. In summary, our models predict an 
increase of the importance of the FeCZ for cooler surface tem- 
perature or lower surface gravity, for higher luminosity, and for 
higher metallicity. 

While in the discussed range of luminosity and effective tem- 
perature, the average convective velocity {vc) is on the order of 
1 to lOkms-' for the FeCZ, we found that the average convec- 
tive velocity {vc) in the He convective zone is always very low 
(< 1 kms-'). Convection due to hydrogen recombination is ab- 
sent; this dominates at lower effective temperatures than the ones 
studied here. 



For our grid of stellar evolution models, we map the average 
convective velocity of the FeCZ (Eq. |6]l in the HR diagram for 
the three different metallicities (see Fig.|9] and Sect. l4.1.2] i. This 
figure displays the three qualitative trends of the iron zone we 
have just described. 

- For given luminosity and metallicity, the average convec- 
tive velocity near the upper boundary of the FeCZs increases 
with decreasing surface temperature. The convection zones 
are located deeper inside the star (in radius, not in mass), 
and the resulting larger pressure scale height leads to higher 
velocities. At solar metallicity and 10^ (i.e. roughly at 
20 Mq) the velocities increase from just a few kms"' at 
the ZAMS to more than lOkms-' in the supergiant regime, 
where {vc) - 2.5 kms-' is achieved at Teff - 30 000 K. At 
the lowest considered metallicity, the FeCZ is absent at the 
ZAMS at 10^ Lq, and a level of {vc) - 2.5 kms-' is only 
reached at Teff - 20 000 K. 

- For fixed effective temperature and metallicity, the iron zone 
convective velocity increases with increasing luminosity, 
since a larger flux demanded to be convectively transported 
requires faster convective motions. Figure |9] in Sect. 14.1.21 
also shows that there are threshold luminosities below which 
FeCZs do not occur, i.e., below about IO-^'^Lq, lO-^'^L©, and 
10'*^ Lq for the Galaxy, LMC and SMC, respectively. 

- The FeCZs become weaker for lower metallicities, since due 
to the lower opacity, more of the flux can be transported by 
radiation. The threshold luminosity for the occurrence of the 
FeCZ quoted above for Z=0.02 is ten times lower than that 
for Z=0.004. And above the threshold, for a given point in 
the HR diagram, the convective velocities are always higher 
for higher metallicity. 

3.1. Rotating models 

We considered two 20 M© models with metallicity Z=0.02, one 
rotating at birth with an equatorial velocity of 250 kms-' (cor- 
responding to about 35% of the critical velocity) and one with 
350 kms-' (about 50% of the critical velocity). The evolution of 
the radial extent of sub-surface convection in the rotating mod- 
els is very similar to the non-rotating case shown in Fig.|2] Also 
the convective velocities inside the FeCZ change only a few per- 
cent between rotating and non-rotating models, even if the rotat- 
ing models show slightly higher convective velocity peaks (see 
Fig. nil. We conclude that rotation is not significantly affecting 
the structure and the properties of sub-surface convection in the 
vast majority of OB stars. 
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Fig. 4. Convective velocity in the FeCZ as function of radial dis- 
tance from the stellar surface. The dotted line corresponds to a 
non-rotating 20 Mq model at Z=0.02, while the dashed and solid 
lines refer to the same model rotating at birth with 250 km s"' 
and 350 kms"' respectively. The values correspond to mod- 
els having the same effective temperature (log T - 4.339) and 
very similar luminosity (log L/Lq = 5.04 for the non-rotating 
model and log L/L© - 5.03 for the rotating ones). The gray 
band shows the upper 1.5 pressure scale heights of the FeCZ, 
which is the region considered for the computation of {vc), cf. 
Eq.|6] Convective velocities in the He convection zone are much 
lower than 1 km s ' and are not visible in this plot. 



As pointed out bv lMaeder et al.l (|2008'), the effects of rotation 
on sub-surface convection become substantial for stars rotating 
close to critical velocity. While stars rotating with such high ve- 
locities exist (e.g. Be stars), their number is modest. The study of 
sub-surface convection in these very fast rotators is interesting, 
but may require 2-dimensional stellar models, which is beyond 
the scope of this paper. 

4. Comparison with observations 

In the following, we investigate the idea that these sub-surface 
convection zones might be related to observable phenomena at 
the stellar surface. In particular, we investigate potential connec- 
tions with microturbulence in massive stars, and discuss whether 
small scale or large scale clumping in massive star winds, mag- 
netic fields, and non-radial pulsations could be related to sub- 
surface convection. For each point, we first briefly discuss the 
theoretical motivation, and then the corresponding observational 
evidence. 

4.1. Microturbulence 

4.1.1. Theoretical considerations 

The convective cells in the upper part of a convection zone excite 
acoustic and gravity waves that propagate outward. The genera- 
tion of sound w aves by turbulent motions was first discuss ed by 
lLighthilll(ll952l) and extended to a strati fied atmosphere bv lSteinI 
1961) and iGoldreich & Kumad(ll990h . In a stratified medium. 



gravity acts as a restoring force and allows the excitation of grav- 
ity waves. For both acoustic and gravity waves, the most impor- 
tant parameter determining the emitted kinetic energy flux is the 
velocity of the convective motions. This is why, in the follow- 



ing, we use the average convective velocity {vc) as the crucial 
parameter determining the efficiency of sub-surface convection. 

Goldreich & Kumar ( 1990) showed that convection excites 
acoustic and gravity waves, resulting in maximum emission for 
those waves with horizontal wave vector kk ~ 1 /Hp c and angular 
frequency oj ~ Uc/Hp_c, where now Uc and Hp ^ are evaluated at 
the top of the convective region. They calculated that the amount 
of convective kinetic energy flux going into acoustic and gravity 
waves is 



and 



(7) 



(8) 



respectively, where we take Fc ~ pd^c)^ and is the Mach 
number in the upper part of the convective region. Since con- 
vection in our models is subsonic, gravity waves are expected 
to extract more energy from the convective region than acoustic 
waves. These gravity waves can then propagate outward, reach 
the surface and induce observable density and velocity fluctua- 
tions (Fig.|5]l. 

The Brunt- Vaisala frequency in the radiative layer above 
the FeCZ is about mHz. Molecular viscosity can only damp 
the highest frequencies, while wavelengths that will be resonant 
with the scale length of the line forming region should not be 
affected (see e.g. Lighthill 1967). This is the case for the grav- 
ity waves stochastically excited by convective motions: they can 
easily propagate through the sub-surface radiative layer, steep- 
ening and becoming dissipative only in the region of line forma- 
tion. 

Again, multi-dimensional hydrodynamic simulations would 
be the best way to compute the energy loss of these waves during 
their propagation through the radiatively stable envelope above 
the FeCZ, but this is beyond what we can presently do. We can, 
however, obtain an upper limit to the expected velocity ampli- 
tudes at the stellar surface, where we only consider the energy 
transport through gravity waves. The kinetic energy per unit vol- 
ume associated with the surface velocity fluctuations must 
be comparable to or lower than the kinetic energy density as- 
sociated with the waves near the sub-surface convection zone. 



Eg ~ McPcivcf, or 



> 1, 



(9) 



where p^ is the density at the top of the convective region and ps 
is the surface density, and is the surface velocity amplitude. In 
this ratio we only consider energy density since the volume of 
the line forming region is comparable to the volume of the upper 
part of the convective zone. Therefore, we expect 



Ps 



(10) 



In our models with well developed FeCZs, -sjMcPc/ps - 1 (or- 
der of magnitude), and thus and {u^} should be on the same 
order of magnitude. It is difficult to estimate the typical corre- 
lation length of the induced velocity field at the stellar surface, 
but a plausible assumption is that it is about one photospheric 
pressure scale height. Hp s, given the proximity of the FeCZ to 
the surface and the fact that the horizontal wave vector of the 
emitted waves is fc/, ~ 1 /Hp ^. 
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Fig. 5. Schematic representation of the physical processes con- 
nected to sub-surface convection. Acoustic and gravity waves 
emitted in the convective zone travel through the radiative layer 
and reach the surface, inducing density and velocity fluctuations. 
In this picture microturbulence and clumping at the base of the 
wind are a consequence of the presence of sub-surface convec- 
tion. Buoyant magnetic flux tubes produced in the convection 
zone could rise to the stellar surface. 



4.1 .2. Observations 

The microturbulent velocity ^ is defined as the microscale non- 
thermal component of the gas velocity in the region of spectral 
line formation: 



AAn = - 



l-RT 



(11) 



Assuming that the gas in this zone has a temperature only 
slightly diff'erent from the effective temperature, one finds em- 
pirically that the observed Doppler widths AA p cannot be ac - 
counted for by the thermal motions alone (e.g. ICowlevlll97(il) . 
Regardless of which physical mechanism causes microturbu- 
lence, the process of spectral line fitting yields values of f in hot 
massive stars between and about 20 km s"' . In contrast, macro- 
turbulence corresponds to velocity fluctuations which are coher- 
ent on a length scale larger than the radial extent of line forming 
regions. If indeed the length scale of the photospheric velocity 
fluctuations induced by the iron convection zone are on the order 
of the photospheric pressure scale height, then this length scale is 
also comparable to the radial extent of line forming regions, and 
it is difficult to decide whether the velocity fluctuations would 
be manifested as micro- or as macroturbulence, or both. Below, 
we compare our model predictions only to the case of microtur- 
bulence since this is the empirical parameter most extensively 
available in the literature. 

Photospheric microturbulence is routinely required, e.g., to 
derive consistent surface abundances for one element from dif- 
ferent photospheric absorption li nes throug h stellar mode l atmo- 
spheres (among many other s .RoUeston et al.J996; Hibbins et al.l 
ll998l:IVrancken et al]|200(]|) . Unfortunately, difi'erences in phys- 
ical assumptions or atomic physics can require somewhat dif- 
ferent microturbulent velocities for the same star in different 
st udies. Here, we restr ict our detailed c omparison to the data 
of iTrundle et"ai] (|2007|) and lHunter etal] (f2()08b,) from the ESQ 
VLT-FLAMES Survey of Massive Stars jEvans et al.l |2005|) . 
since it comprises the largest available uniformly analyzed data 
set. In Fig|6] we plot the microturbu lent velocities derived for the 
LMC early B type stars analvzed bv lHunter et al.l ( l200 8b) versus 
their projected rotational velocity. The error bar on the derived 
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Fig. 6. Projected rotational velocity usin/ versus photospheric 
microturbule nt velocity ( f or the early B-type stars in the LMC 
analyzed by Hunter et al.l ([2008b). Different symbols refer to 
different luminosity intervals, as explained in the legend. The 
microturbulent velocities ^ have typical uncertainties of about 
+5kms"'. An uncertainty of 10% or ±10kms"', whichever is 
the larger, should be considered for the rotational velocity mea- 
surements. 



microturbulent velocities is usually quite big, ±5kms"', and is 
often comparable to the measured quantity itself. There seems to 
be no positive correlation between ^ and the apparent projected 
rotational velocity v sin /. Though v sin / is plotted and not v itself, 
the lack of a correlation in such a large data set (justifying the 
assumption of random orientation of the sample) argues against 
rotation as an important effect in triggering microturbulence in 
hot stars. To compare microturbulent velocities to properties of 
sub-photospheric convection we use only data obtained for slow 
rotators (i.e. ysin/ < 80kms"') as microturbulent velocities are 
more difficult to measure for faster rotators. 

In Fi g IT] we show t he mi croturbulent velocities for the LMC 
stars of 



Hunter et al.l t2008b) versus the stellar surface grav- 



ity. Trends of the microturbulent veloci ty with log ^ have been 
previously reporte d for hot stars (e.g iGies & LambertI 119921; 
iHunteret all 12007b . The fi gure shows that indeed, for log^ < 
3.2, there is a clear trend. However, the luminosity coding 
in Fig|2] suggests that this trend may be largely produced by 
the increase in convective velocity with increasing luminosity 
(Sect. O. This figure displays a detection threshold of about 
lOkms"' for the microturbulent velocities so in the following 
we restrict the comparison to ^ > 10 km s"' . 

In order to compare these observations to our model pre- 
dictions, we evaluated the ratio of the kinetic energy in the 
form of gravity waves at the surface of the FeCZ to the ki- 
netic energy of the surface velocity field, Eg/Es (Eq. |9l), as- 
suming i)s = lOkms ', in the HR diagram. Fig. [8] shows two 
different iso-contours of this ratio; the stars of the LMC sam- 
ple shown in Fig. |7] are over plotted. Notably, all but one of 
the LMC stars of Fig. |8]with ^ > lOkms ' are found in that 
part of the HR diagram where it is energetically possible that 
the FeCZ-induced gravity waves trigger a significant surface ve- 
locity field (iij > lOkms '). Thus, a physical connection of 
the FeCZ with the observed microturbulent velocities appears 
energetically possible. Moreover, that the iso-contour line of 
Eg/Es = 1 in Fig. |8]almost perfectly divides the observed sam- 
ple in stars with significant (^ > lOkms ') and insignificant 
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Fig. 7. Logarithm of surface gravity versus micro turbulent ve- 
locity ^ for the LMC early B-type stars studied by iHunter et al.l 
(|2008b); only stars with vsini < 80kms ' are considered here. 
Different symbols refer to different luminosity intervals, as ex- 
plained in the legend. The microturbulent velocities ^ have typi- 
cal uncertainties of about ±5 km s"' . For the surface gravity mea- 
surements an uncertainty of ±0. 1 should be considered. 



< 10 km s ') microturbulence is a further indication of such a 
physical connection. 

Figure |9] shows the iso-contours in the HR diagram of the 
average convective velocity from our models in the upper layers 
of the iron convective zone, (vc) (cf.. Sect. 14.1.1b . at the three 
considered metallicities. We have over plotted the microturbu- 
lent vel ocities derived bv lTrundle et alj (l2007l) and lHunter et alj 
(12008b') as filled circles. A gain, we distinguish between sam- 
ple stars with significant > lOkms '; Group A) and in- 
significant < lOkms '; Group B) microturbulent veloci- 
ties. Comparing the plot for the LMC in Fig. |9] with Fig. [8] 
identifies (uc) - 2.5 km s"' as a critical convection velocity to 
be able to trigger microturbulence. Interestingly, the contour of 
(Vc) = 2.5 km s"' in our stellar models forms an almost perfect 
dividing Une between Groups A and B for all three considered 
metallicities. 

In fact. Fig. |9] provides evidence for all three trends found 
in the average convection velocity as function of stellar param- 
eters (cf.. Sect. [3]) to be present also in the empirical data on 
microturbulent velocities. The LMC data shows that in the lumi- 
nosity range 4.5 < log L/Lq < 5.5 microturbulence is found 
only for Teff < 25 000 K. The data for all three metallicities 
clearly suggests a key role of the luminosity, as the stars with 
^ > lOkms"' are the most luminous ones in each sub sample. 
And finally, the stars with high microturbulent velocities are all 
comfortably above our theoretical contour line corresponding to 
(uc) - 2.5 km s"'. As the latter trends toward the upper right 
corner of the HR diagram for lower metallicity, the metallicity 
dependence is also confirmed by the empirical data. 

iLvubimkov et alj (^004) studied microturbulence in a sam- 
ple of 100 Galactic early B stars. Interestingly, they found sig- 
nificant microturbulent velocities (i.e., clearly above 5kms"') 
in the mass range 7...1 1 M© for stars with a relative age on the 
main sequence of J/tms > 0.8, and in the range 12...19M0 for 
f/i"MS > 0.25, but only insignificant microturbulent velocities 
for younger or less massive stars. Again, these results appear to 
agree with Fig.|9]up to a remarkable quantitative level. 
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Fig. 8. Values of the ratio Eg/Es of the kinetic energy in the form 
of gravity waves above the iron convection zone, to the kinetic 
energy of the surface velocity field, as a function of the loca- 
tion in the HR diagram (see color scale). This plot is based on 
evolutionary models between 5 Mq and lOOM© for LMC metal- 
licity. We estimated the ratio Eg/E^ as in Eq. |9] using a value 
Us = lOkms ' for the surface velocity amplitude. Over-plotted 
as filled circles are stars which have photos pheric micro t urbulen t 
velocities ^ derived in a consistent way by Hunter et al.l (l2008bl) . 
Here, we use only data for stars with an apparent rotational ve- 
locity of vsini < 80kms The uncertainty in the determina- 
tion of ^ is typically +5kms"\ which justifies our choice of 
Vs = lOkms"^ Solid white lines are reference evolutionary 
tracks. The full drawn black line corresponds to the zero age 
main sequence. 



In summary, our comparison provides evidence for a physi- 
cal connection of microturbulence in hot star photospheres with 
the existence and strength of a sub-photospheric FeCZ. 

If microturbulence has a physical origin and is not just a 
fudge factor, the pressure and energy terms associated with such 
a velocity field should be includ ed in the calcu l ation s of atmo- 
spheric models of massive stars. iHubenv et al"] (1199 1) have in- 
vestigated part of these effects by accounting for a constant mi- 
croturbulent velocity in the pressure term only. They find that 
for stars with conspicuous ^ values (of 25kms"') the inclusion 
of the pressure term leads to higher values of the surface gravity, 
which can reduce the mass discrepancy for O stars and O-type 
central sta rs of planetary nebula. A similar approach was also 
studied by Smith & Howarth ( 1998). The impact on gravity dis- 
cussed by Hubeny et al. (1991) is likely an upper limit to the 
effect as, first, the f values are in most cases less than 25 km s ' , 
and, second, a positive gradient in the atmospheric ^(r) would 
decrease the pressure gradient due to microturbulence but, to 
date, the radial stratification of the microturbulent velocity in 
the atmospheres of hot massive stars has not been studied in de- 
tail. From a theoretical perspective, investigating ^(r) requires 
hydrodynamic simulations of the stellar atmosphere, including 
the presence of sub-surface convection. 

The mass discrepancy in massive stars is a we l l doc- 
umented problem ( see for example iHerrero et al.l 119921; 



Lennon et al.l 120031: iTrundle & LennonI 120051; iMassev et al] 



i2005[ iMokiem et all 12007"). It is typically found that the dif- 
ference between spectroscopic mass and evolutionary mass is 
most pronounced in supergia nts. In main s e quence stars it may 
not be present at all, but see iHunter et al.l (l2008bl) . Given that 
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Fig. 10. Values of the mass discrepancy (evolutionary mass di- 
vided by spectroscopic mass) as function of microturbulent ve- 
locity in the sample of B stars analysed bv iTrundle et alj (|2007|) 
and Hunter et al. (2008b). Here, we use only data for stars with 
an apparent rotational velocity of v sin / < 80 km s"^ 



microturbulent velocities are highest in supergiants (see Fig. [T]) 
an empirical correlation between mass discrepancy and micro- 
turbulent velocity is to be expected a nd is show n in Fig. [TOl 
using data analysed by iTrundle et alj (l2007h and iHunter et al.l 
(2008b). If indeed microturbulence is related to subsurface con- 
vection and supergiants have intrinsically higher microturbulent 
velocities than dwarfs (see Section[3]) potentially part of the gra- 
dient in Fig. [TOl ma y be explained by the effect discussed by 



4-0 iHubenv et al.1 11991 . 
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Fig. 9. Average convective velocity within 1.5 pressure scale 
heights of the upper border of the iron convection zone in our 
models, as function of the location in the HR diagram (see 
color scale), based on evolutionary models between 5 Mq and 
100 M0 (white lines), for three metallicities corresponding to the 
Galaxy (top panel), the LMC (middle), and the SMC (bottom). 
The full drawn black line corresponds to the zero age main se- 
quence. Over-plotted as filled circles are photospheric microtur- 
bulent velocities ^ derived in a consistent way for hot massive 
stars by Trundle et al. (2007) and Hunter et al.1 (|20b8b). Here, 
we use only data for stars with an apparent rotational velocity of 
ysin; < 80kms The uncertainty in the determination of ^ is 
typically ±5kms '. 



4.2. Non-radial pulsation 

4.2.1. Theoretical considerations 

In our discussion thus far we have considered only the prop- 
agation of running waves, it is possible that the stochastic 
convective motions can also excite standing waves, i.e. high- 
order non-radial pulsations. For example, stochasti c excitation 
is thought to be the cause of the Solar oscillations ( UMchll 19701: 
Leibacher & Stein 1971). It may thus be possible that the FeCZ 
excites non-radial pulsations in hot early-type stars. 

Several classes of OB star models are found to be linearly 
unstable against non-radial pulsations, among which are the 
13 Cephei stars and the slowly pulsating B stars (e.g., Dupretl 
'200 1 1 |Pamyatnyklv 1999). The key ingredient required for the 
pulsational instability is the iron opacity peak described in 
Sect. 12] As convection is not required to produce the pulsations 
in thes e models, it is not co nsidered in detail as excitation mech- 
anism (Dziembowski 2008). It is conceivable that the convective 
excitation could modify the predicted pulsation spectrum and/or 
extend the instability region of certain linear instabilities. The 
convective kinetic energy flow into waves could be predomi- 
nantly directed into those modes for which instability is pre- 
dicted in the models. In certain parts of the HR diagram, one 
may thus suspect an intricate connection between the occurrence 
of a sub-photospheric iron convection zone and the properties of 
non-radial pulsations. 

Non-radial pulsations have also been considered as the ori- 
gin of various observed small scale (e.g., line profile variabil- 
ity, iFuUerton et al.ll996llI997l) and large scale phenomena (e.g.. 
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so ca l led discrete absorption comporients iPrinia & HowarthI 
fT988t iMassa et all fT995L iKaoer et all \l99% iPrinia et alJ 120021) 
at the surface or in the wind of massive OB stars. Non-radial 
g-mode pulsations were also recently proposed as t he origin of 
obser vable macroturbulence in massive B type stars jAerts et al.l 
I2008h . In Fig. [m we compare the regions where strange mode, 
g-mode, and p-mode pulsations are predicted to occur in the 
HR diagram with the region where our models predict a strong 
FeCZ. Pulsations appear to be almost ubiquitous when all types 
of variables are accounted for The strange mode pulsators are 
predicted to cover the HR diagram at high luminosity, where 
w e plotted only the pre dictions for the radial strange modes 
of lKiriakidis et alJ (Il993h : high-order non-radial strange modes 
seem to be omnipresent as well for stars above 40 M© or so 
jGlatzel & Mehren"l996'). Non-radial g-mode pulsators are pre- 
dicted by Saio et al. (2006) in the B supergiant region. And 
radial and low order non-r adial p-m odes are pred icted for the 
P Ceph ei r egime bv .Deng & Xiond (12001 ) and by Pamvatnvkhl 
(119991) and lSaio et all (120061) ' for a considerably larger region in 
the HR diagram. At lower metallicity, many of the predicted ar- 
eas in the HR diagra m are smaller (cf., Kiriakidis et al. 1993; 
iDeng & Xiongll200ll) but the general picture is still incomplete. 

4.2.2. Observations 

Observationally, the classical jS Cephei sta rs are concentrated 
in the region predicted by Deng & Xiong dStankov & Handleil 
I2005h . while the B supergia nt non-rad ial g-mode pulsators 
overlap with the prediction of ISaio et alj (2006) but extend t o 
an effective temperature of ~ 10 000 K CLefever et al.1 12007[) . 
Pulsations are al so found for the most lumi nous stars (e.g., the a 
Cvgni-variables: Ivan Leeuwen et al.lll998l) . but there is now no 
clear evidence for strange mode pulsators. Comparing the pre- 
diction for the FeCZ with that for pulsational instability (FigfTTTl 
shows two things. Firstly, the FeCZ-region is much larger than 
any region for a particular pulsational instability. Thus, distin- 
guishing whether a certain observational feature is caused by a 
particular pulsational instability by the FeCZ might, in principle, 
be possible, since the area in the HR diagram where the latter 
occurs but the pulsational instability does not is relatively large. 
Secondly, some regions exist where (so far) no pulsations are 
predicted but the FeCZ in our models is strong, or where, vice 
versa, pulsations are predicted but the FeCZ is weak or absent. 

Comparing Fig. [TT] with Fig.|9] where we show the observa- 
tions of microturbulence and the FeCZ predictions, it is unlikely 
that microturbulence is associated with a particular pulsational 
instability. Strong microturbulence is observed at too low a lumi- 
nosity to be attributable to strange mode pulsations alone, while 
p-mode pulsators are found where microturbulence seems not 
to occur Concerning the g-mode pulsators the situation is less 
clear Fig. [TTI shows that, at solar metallicity, g-mode pulsations 
for post-main sequence stars are expected only in a rather narrow 
luminosity interval. Unfortunately, the five Galactic stars shown 
in Fig. 5 for which strong microturbulence is derived are all in- 
side this luminosity range, so they cannot distinguish between a 
pulsational or FeCZ origin of microturbulence. However, look- 
ing at the LMC data, stars above the g-mode luminosity up- 
per limit with microturbulence are found; whether or not cor- 
responding stellar models are g-mode unstable is currently not 
known. A connection of microturbulence with non-radial pulsa- 
tions is thus not impossible, but it is also not very likely. 

Comparing Fig. [TT] with the discrete absorption compo- 
nents (DACs) found in 200 Galactic O stars above ~ 20 by 
iHowarth & Prinial d 19890 all the way to the zero age main se- 
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Fig. 11. The plot shows regions of the HR diagram where pul- 
sational instabilities are predicted, compared to our calcula- 
tions for the occurrence of iron convection. The cloudy region 
marks the presence of iron sub-surface convection with {vc) > 
2.5 kms while the dotted, blue line divides regions of the HR 
diagram where iron convection is present (above) from regions 
where it is absent (below). Different modes of instabilities are 
shown with different colors and different contour line styles, as 
explained in the legend. Evolutionary tracks between 7 Mq and 
40 Mq are plotted as a reference. The straight, full drawn black 
line corresponds to the zero age main sequence. The Humphrey- 
Davidson limit is also plotted for reference (top-right corner). 



quence, seems to argue against non-radial pulsations as the ori- 
gin of the DACs phenomenon (see also Sect. 14.4b . 

4.3. Wind dumping 

4.3.1. Theoretical considerations 

Observational evidence exists for stellar wind inhomogeneities 
on small and on large scales. While the latter will be discussed 
in Sect. 14.41 here we consider only small scale wind structure, or 
wind clumping. In Sect. 14. 11 we discussed that waves produced 
by the FeCZ could lead to velocity fluctuations at the stellar sur- 
face. In order to induce wind clumping, those waves should in- 
duce density fluctuations at the stellar surface. Through the oc- 
currence of porosity or shifts in the ionisation balance of the 
gas the mass-loss rate may be affected. For this to happen, the 
amplitude of the velocity fluctuations at the surface should be 
on the same order of the sound speed. Alternatively, the velocity 
fluctuations might directly affect the local mass-loss rate through 
the Doppler effect, if the amplitude of the velocity fluctuations 
is on the same order as the speed of the wind flow, which, at 
the base of the wind, is approximately the sound speed. As the 
sound speed at the surface in our massive main sequence models 
is on the order of a few times 10km s ', we consider here those 
stellar models potentially capable to produce wind clumping for 
which the convective velocities in the upper part of the FeCZ 
{vc) > 2.5 kms"', as this allows energetically to have surface 
velocity amplitudes above ~ 10 kms"' (cf. Sect. 14. II ). 

Assuming the horizontal extent of the clumps to be compa- 
rable to the sub-photospheric pressure scale height H^, we may 
estimate the number of convective cells by dividing the stellar 
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surface area by the surface area of a convective cell finding that 
it scales with (R/Hp)^. For our main sequence O star models in 
the mass range 20-60 Mq, we find pressure scale heights in the 
range 0.04-0.24 R©, corresponding to a total number of clumps 
in the range 6 xlO^ - 6 x 10**. In principle, this might be testable 
through linear polarization variability measuremen ts, which can 
probe wind asphericity at the very base of the wind dPavies et alj 
I2007h . 

4.3.2. Observations 

Evidence has been accumulating that the winds of massive stars 
may be subject to small scale clumping. So far this is best doc- 
umented for Wolf-Rayet (WR) stars, where line variability on 
time scales of minutes to hours is thought to constitute direct ev- 
idence of outflows that are clum ped already in the accele ration 
zone near the base of the wind dLepine & Moffatlll999() . This 
clumping may be part of the explanation for t he wealth of in- 
tricat e detail see n in nebulae around WR stars jGrosdidier et alj 
11998 ). Recently. iLepine & Moffa t' (2008) reported spectr oscopic 
variab ility in the Of supergiants ^Pup (see also Eversbe rg et alj 
119981) and HD 93 129 A. The amplitude of the variation (at the 
1-3% level) is similar as in WR stars supporting the notion that 
clumping is not restricted to WR stars. 

Indeed, evide nce that O star winds are clumped is given by, 
among others. Puis et alj ( |2006.) . These authors investigate the 
clumping behavior of the inner wind (inside about two stellar 
radii) relative to the clumping in the outer wind (beyond tens 
of stellar radii) of a large sample of supergiant and giant stars. 
They find that in stars that have strong winds, the inner wind 
is more strongly clumped than the outer wind, whereas those 
having weak winds have similar clumping properties in the inner 
and outer regions. Their analysis only allows for such a relative 
statement. In principle, for weak winds the outer part could be 
homogeneous. If so, weak winds are not clumped. In any case, 
strong winds - identified as such if Ha is seen in emission - are 
clumped near the base of the wind. A measure of the degree 
of clumping is the clumping factor /d - (p^)/(p)^ > 1 where 
angle brackets denote (temporal) average values (e.g. lPuls et alj 
120061) . 

Apparently, this type of radial behavior is not consistent 
with hydrodynamical predictions of the intrinsic, self- excited 
line-driven instability iRunacres & Owockil llOOl l2005h . Such 
models predict a lower clumping in the inner wind than the 
outer wind. Moreover, if there was any dependence on wind 
density predicted at all, optically thin winds should be more 
strongly clumped than optically thick winds dOwocki & Pulsl 
11999 ^: Puis et al. 2006). Therefore, the findings on the radial 
clumping behavior in O stars may point to an additional exci- 
tation mechanism of wind structure. 

J'ig. [12] shows that the O stars investigated by iPuls et afl 
populate the regime in the HR diagram in which our 
models predict the average convective velocity in the top part 
of the FeCZ to change from a few to over 2.5 kms ', indicat- 
ing that surface velocity fluctuations on the order of the lo- 
cal sound speed are possible (cf. Sect. 14.3. It . Though the part 
of the HR diagram that is covered by the sample is limited 
(4.46 <> log Teff <: 4.66; 5.29 <: log L/Lq $ 6.26), the trend 
is such that stars with relatively strong clumping in the inner 
winds are in a regime where (ue) is higher. A correlation be- 
tween clumping at the base of the wind and (uc), i.e., between 
wind clumping and the properties of the FeCZ, appears there- 
fore possible. To further test the idea that the FeCZ produces 
wind clumping at the wind base for sufficiently luminous and 
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Fig. 12. Average convective velocity within 1.5 pressure scale 
heights of the upper border of the iron convection zone in our 
models, as function of the location in the HR diagram (see 
color scale), based on evolutionary models between 5 M© and 
120Mo (white lines) at solar metallicity. The full drawn black 
line corresponds to the zero age main sequence. Over-plotted as 
filled circles are observations of the clumping factor f " (see text 
for de finition) in the winds of O stars, according to [Puis et alj 
(12006!). The data shown here corresponds to objects with well- 
constrained clumping parameters. Note the different luminosity 
range with respect to Fig.|9l 



cool stars it would be desirable to derive the radial clumping pro- 
files for cooler (i.e. B-type) stars. If correct, such stars, both the 
ones with weak and strong winds, should have relatively strong 
clumping at the base of the wind. 

To derive the spatial scale of the wind clumps from linear 
polarimetry has not yet been possible for main sequence OB 
stars. A limitation is that this technique requir es very high 
signal-to-noise observations (see discussion in iHarries et alj 
2002). Luminous Blue Variables (LBVs) however provide a 
more appropriate category of objects to test wind clump sizes, 
because of a combination of high er mass-loss rates, a nd lower 
wind velocit i es tha n for O stars dPavies et alj l2005l) . Indeed, 
iDavies eta D (200"/) show that in order to produce the observed 
polarization variability of P Cygni, the wind should consist 
of about ~ 1000 clumps per wind flow- time (t// = Ri,/uoa)- 
To see whether this observational result is compatible with 
sub-surface convection causing wind clumping, we considered 
the sub-surface convective regions of a massive star model with 
global properties similar to those of P Cygni (initial mass 60 
Mo, log(L/L0 = 5.9, and log T = 18 000 K ). As a result of 
the lower gravity, the pressure scale height in the FeCZ in this 
model is about 4 R©, which is much bigger than in our O star 
models. Consequently, the same estimate for the number of 
clumps as done for the main sequence models in Sect. 14.3.1! 
yields about 500 clumps per wind flow time, a number which 
is quite comparable to that derived for P Cygni observationally 
(about lO-' clumps per wind flow time). 



Finally, iFullerton et alj (Il996l) have conducted a spectro- 
scopic survey of O stars and observed intrinsic absorption line 
profile variability (LPVs) for about 77% of their sample. They 
report an increase of incidence and amplitude of variability with 
increasing stellar radius and luminosity, as well as no statistically 
significant line profile variability for dwarfs earlier than 07. 
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While Fullerton et al. attempt to relate their fin dings to the pre- 
diction s of strange-mode pulsation in O stars bv lKiriakidis et alj 
(Il993h . a comparison of their results (see their Fig. 13) with the 
occurrence of sub-surface convection as depicted in Fig.|9]indi- 
cates the possibility of a physical connection between line profile 
variability and sub-surface convection in O stars. 

4.4. Magnetic fields 

AAA. Fields from iron convection zones 

In solar-type stars, surface convection zones modified by the 
stellar rotation are thought of being capable of producing a 
magnetic field through the so called gQ-dyn amo (iParkerll 19751; 
ISpiegel & WeisslfigSOHSpTegel & ZahnlfT992h . The FeCZ in our 
massive main sequence stellar models has a spatial extent sim- 
ilar to the Solar convection zone, although its mass is much 
smaller, and OB stars are rapid rotators, so it is possible that 
a dynamo may also work in the envelopes of OB stars. If so, 
the magnetic field may be able to penetrate the radiatively sta- 
ble layer above the FeCZ, and dynamically significant field 
str engths might be achievab l e. To this end, we follow the model 
bv lMacGregor & Cassinelhl (l2003i) for the rise of buoyant mag- 
netic flux tubes generated at the edge of the convective core 
of massive stars through their radiative envelope and apply this 
model to the FeCZ and the overlying radiative layer The mag- 
netic field strength Bq in the iron convection zone is estimated 
assuming equipartition of kinetic energy density and magnetic 
energy density inside the convective layers: 

Bq ^ 2vc yfnp, (12) 

which, for our 60 Mq star at Z=0.02, reaches about Bq - 2700 G 
inside the iron convective zone. The surface field is then 
obtained by multiplying this number with the ratio of the sur- 
face density ps and the density in the FeCZ po, i.e. B^ ^ 
BoPs/po - 60 G. Similarly, for the 20 Mq model at Z=0.02 we 
obtain Bo - 1400 G and 5, ^ 10 G. Although at the surface, 
the magnetic pressure in the flux tubes is only on the order of 
a few percent of the total pressure, it is on the same order as 
the gas pressure and could thus lead to considerable horizon- 
tal density diff'erences. Compar ed to the situation envisioned by 
iMacGreg or & Cas sinellil (l2003h . who found that the rise time of 
the flux tubes from the edge of the convective core to the stellar 
surface can be comparable to the ma in sequence life time (but 
see also MacDonald & Mullanll2004l) . the rise time of the flux 
tubes from the FeCZ to the surface is much shorter And while 
the initial magnetic field strength at the edge of the convective 
core can be considerably higher than our values of Bq, the sur- 
face fields obtainable from the sub-surface convection zones are 
higher, due to the much lower density contrast between convec- 
tion zone and surface in this case. 

As a consequence, even though we are far from a detailed 
picture, it seems conceivable that the FeCZs in massive main se- 
quence stars produce localized magnetic fields at their surface. 
The interaction of the stellar wind with the localized surface 
magnetic fields could enhance the rate at which the wind induces 
a loss of stellar angular momentum. Furthermore, co-rotating 
density patterns in the outflowing wind could be produced by 
these local magnetic spots. 

Rotation may play an important role in the dynamo process, 
possibly resulting in the appearance of stronger fields at the sur- 
face for faster rotating stars. To estimate this effect, a dynamo 
model accounting for the differential rotation needs to be im- 



plemented in the stellar evolution calculations. This will be the 
discussed in a subsequent paper. 

4.4.2. Observations 

Surface magnetic fields have been linked to several observed 
phenomena in OB stars, e.g. discret e abso rption compo nents 
(DACs) in U V resonance fine s (e.g . . iPrinia & Howarth 119881: 
iMassa et alJ[T995l: iKaper e"tanil997l [Prinia et al..,2002) . which 
are thought to diagnose la r ge scale coherent wind anisotropics 
dCranmer & Owockil ll996': Lobe l & Blommd l2008'), or t he less 
coherent fin e profile variability mentioned above ( Fullerto n et al.l 
[19961 fT997h . Also non-thermal X-ray emission of OB main se- 
quence stars has been proposed to relat e to surface magnetic 
fields (e.g., iBabel & Montmerld 119971: lud-Doula & Owockil 
llOOl) . 

A connection of the FeCZ in massive stars with the phe- 
nomena mentioned above has not yet been considered. However, 
such a connection becomes testable through our results. While 
in our comparison to observed microturbulence presented above 
we discussed when sub-surface convection may lead to de- 
tectable surface velocity fluctuations, the presence of surface 
magnetic fields may simply depend on whether an FeCZ is 
present in the star or not. Looking at Fig. |9] we see that in 
our models the FeCZ is absent for luminosities below about 
10^ 2 10^ ') Lq, and 10^* ^ Lq for the Galaxy, LMC and SMC, 
respectively. If DACs or line profile variability were produced 
by magnetic flux tubes generated in the FeCZ, those phenomena 
would not be expected f or OB stars below those luminosities. 
iHowardi & Prinial (Il989l) find DACs in nearly all O stars (97%) 
in a large Galactic sample, with log L/Lq > 4.5 and with ef- 
fective temperatures as high as the zero-age main sequence val- 
ues of stars above ~ 20 M©. Since those stars are well above 
the luminosity threshold for the occurrence of the iron convec- 
tion at Galactic metallicity, these observations do not exclude 
DACs being due to FeCZ induced B-fields. Also, all eleven 
early B supergiants with DACs in the sample of Prinia et alj 
(2002) are predicted to have strong FeCZ by our results. Notably, 
between about 20 M© and 40 M©, stars close to the zero-age 
main sequence are not predicted to be pulsationally unstable (cf. 
Fig. fTTT i. which may be in conflict with pulsations as the origin 
for DACs. 



4.4.3. Otiier types of fields 

It may be interesting to briefly compare the expectation from 
surface magnetic fields produced via the FeCZ to that for fields 
produ ced by other means. Surface fields produced by convective 
cores (ISchuessler & Paehler"l978'; 'Charb onneau & MacGregod 
2001; MacGregor & Cassinelli 2003) have been proposed to re- 
late to the same phenomena as those mentioned above, even 
if for massive stars the buoyant rise of magnetic fields from 
the convective core seems to be unhkely dMacDonald & MuUanl 
2004). In contrast to the sub-surface FeCZ, convective cores are 
prevalent in all stars above about 1.2 M©. It has been found that 
the longer lifetime of stars of lower mass may favor the drift of 
fields produced in the core to the surf ace (iSchuessler & Paehlerl 
119781; iMacGregor & Cassinellil l2003h . Therefore, the expected 
trend is opposite to that found for fields produced by the FeCZ, 
where surface fields may occur only for stars above a critical 
mass (or luminosity), and stronger fields are found for more mas- 
sive stars. 



12 



M. Cantiello et al.: Sub-surface convection in hot stars 



On the other hand, in contrast to fields from the FeCZ, 
magnetic flux tubes produced in the core may carry CNO- 
processed material to the surface. This might thus constitute a 
mechanism to ex plaining nitrogen enrich r nent in slowly rotat - 
ing early B stars (iMorel et al.ll2006l l2008t iHunter et al.ll2008al) . 
Strong fossil magnetic fiel ds are thought to persist in only 
a fraction of massive stars dPerrario & W ickramas inghg 120051: 
iBraithwaite & Spruit 2004), and may lead to, among other phe- 
nomena, highly anomalous surface chemical compos itions, wind 
confinement, and var iable X-ray emission (e.g., IWade et aTl 
120061; iTownsend et an i2005). Those strong features can clearly 
not be produced by fields originating from the FeCZs. 

Finally, magnetic fields produced in differentially rotat- 
ing massive stars by the Spruit- Taylor dynamo dSpruitl l2b02h 
may tran sport angular momentum and chemical species (cf., 
iHeger et al. 2005). These fields are predominantly of toroidal ge- 
ometry and would quickly decay near the stellar surface, and are 
thus not tho ught to lead to observab l e field s at the stellar surface 
(but see also lMuUan & MacDonaldl d2005l) ). 



5. Concluding remarks 

Hot luminous stars show a variety of phenomena at their photo- 
sphere and in their winds which still lack a clear physical in- 
terpretation at this time. Among these phenomena are photo- 
spheric turbulence, spectral line variability (DACs and LPVs; 
see Sect. |4|i, wind clumping, variable or constant non-thermal 
X-ray and radio emission, chemical composition anomalies, and 
intrinsic slow rotation. In the previous section, we argued that the 
kon convection zone could be responsible for various of these 
phenomena. 

We proposed in Sect. l4.1l that a physical connection may ex- 
ist between microturbulence in hot star atmospheres and a sub- 
surface FeCZ. The strength of the FeCZ is predicted to increase 
with increasing metallicity Z, decreasing effective temperature 
T and increasing luminosity L (Sect. jSjl, and all three predicted 
trends are reflected in the observational data. This suggests that 
microturbulence corresponds to a physical motion of the gas in 
hot star atmospheres. This motion may then be connected to 
wind clumping (Sect. 14. 3t , since the empirical microturbulent 
velocities are comparable to the local sound speed at the stellar 
surface. In order to verify such a picture, multi-dimensional cal- 
culations of the FeCZ and the radiative layers above, including 
the stellar atmosphere, are required — similar to the recent gen- 
eration of atmosphere models for cool stars (e.g.. lAsplund et akl 
fT999tlWedemever et al.ll2004l) . 

In Sect. 14.41 we proposed that the FeCZ in hot stars might 
also produce localized surface magnetic fields, in Galactic stars 
for luminosities above ~ 10^'^ Lq. This could explain the occur- 
rence of DACs (discrete absorption components in UV absorp- 
tion lines), also in very hot main sequence stars for which pulsa- 
tional instabilities are not predicted. We further argued that there 
may be regions of the upper HR diagram for which the presence 
of the FeCZ influences, or even excites, non-radial stellar pulsa- 
tions (Sect.|4r2]i. 

The FeCZ could also turn out to directly affect the evolution 
of hot massive stars. If it induces wind clumping, it may alter 
the stellar wind mass-loss rate. Such a change would also in- 
fluence the angular momentum loss. In addition magnetic fields 
produced by the iron convection zone could lead to an enhanced 
rate of angular momentum loss. These effects become weaker 
for lower metallicity, where the FeCZ is less prominent or ab- 
sent (see Sect. [3]). 



Finally, we note that the consequences of the FeCZ might 
be strongest in Wolf-Rayet stars. These stars are so hot that the 
iron opacity peak, and therefore the FeCZ, can be directly at the 
stellar surface, or — to be more precise — at the sonic point 
of the wind flow dHeger&Langeilll996l) . This may relate to 
the ve ry strong clumping found observationally in Wolf-Rayet 
winds (L epine& Moffat ' 1999*: 'M archenko et al. F2006). and may 
be required for an unde rstanding of the very high mass-loss 
rates of Wolf-Rayet sta rs [ Eichler et al.ll 995l:lKato & Ibenll992l: 
lHeger&Langeil(l996l) . 
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Appendix A: Models 

Table A.l. Outermost 4 Rq of a 60 Mq model at solar metallicity (Z=0.02). 
The table shows physical variables at t=2.37xl0* years. Columns contain the 
progressive grid point of the model, the status (R=radiative, C=convective), the 
optical depth t, the opacity x, the density p, the radius R, the value M, - 
(where M, is the total stellar mass and Mr is the mass coordinate), the tempera- 
ture T, the convective velocity Vc and the local sound speed Cj. All the values are 
in cgs units if not otherwise specified. 
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1100 
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6.4257E-09 


20.5502 


5.98E-06 


1.504E-I-05 


62.33 


172.73 


1101 


c 


3.82E+02 


1.14352 


6.5118E-09 


20.5989 


5.70E-06 


1.483E+05 


59.95 


167.44 


1102 


c 


3.61E+02 


1.12949 


6.6089E-09 


20.6392 


5.46E-06 


1.464E+05 


57.33 


162.39 


1103 


c 


3.41E-H02 
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6.7225E-09 
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5.23E-06 


1.444EH-05 


54.47 


157.26 


1104 


c 


3.21E-I-02 
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5.00E-06 


1.424EH-05 


51.18 


152.06 


1105 


c 


3.02E-I-02 


1.08331 


7.0035E-09 


20.7520 


4.77E-06 


1.404E-I-05 


47.43 


146.69 


1106 


c 


2.83E-H02 


1.06508 


7.1717E-09 


20.7873 


4.55E-06 


1.382E+05 


43.00 


141.22 


1107 


c 


2.65E-h02 


1.04501 


7.3455E-09 


20.8216 


4.32E-06 


1.360E+05 


37.81 


135.85 


1108 


c 


2.48E+02 


1.02372 


7.5117E-09 


20.8535 


4.10E-06 


1.338E+05 


32.08 


130.70 


1109 


c 


2.31E-I-02 


1.00600 


7.6295E-09 


20.8846 


3.88E-06 


1.319E-H05 


26.60 


126.77 


1110 


c 


2.22E-I-02 


0.99256 


7.7030E-09 


20.9025 


3.76E-06 


1.306E-I-05 


22.40 


123.99 


1111 


c 


2.13E+02 


0.97901 


7.7615E-09 


20.9201 


3.63E-06 


1.292E+05 


19.15 


121.31 


1112 


c 


2.03E-H02 


0.96557 


7.8026E-09 


20.9376 


3.50E-06 


1.277EH-05 


16.16 


118.76 


1113 


c 


1.94E-I-02 


0.95303 


7.8243E-09 


20.9549 


3.38E-06 


1.263E-I-05 


13.55 


116.42 
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Grid point 


STAT 


T 


X 


P 


/?[Ro] 


M. - M, [g] 


T[K] 


Vc [km s"^] 


Cs [km s ^] 


1114 


c 


1.86E+02 


0.94136 


7.8283E-09 


20.9708 


3.26E-06 


1.250E+05 


11.38 


114.27 


1115 


c 


1.78E+02 


0.92946 


7.8150E-09 


20.9866 


3.15E-06 


1.236E+05 


9.54 


112.21 


1116 


c 


1.70E+02 


0.91738 


7.7833E-09 


21.0024 


3.03E-06 


1.222E+05 


7.92 


110.24 


1117 


c 


1.63E+02 


0.90353 


7.7240E-09 


21.0183 


2.91E-06 


1.206E+05 


6.47 


108.10 


1118 


c 


1.53E+02 


0.88804 


7.6268E-09 


21.0386 


2.77E-06 


1.187E+05 


5.08 


105.83 


1119 


c 


1.44E+02 


0.87366 


7.5046E-09 


21.0592 


2.62E-06 


1.169E+05 


3.90 


103.80 


1120 


c 


1.36E+02 


0.86054 


7.3626E-09 


21.0776 


2.49E-06 


1.152E+05 


3.00 


102.00 


1121 


c 


1.28E+02 


0.84800 


7.1942E-09 


21.0964 


2.36E-06 


1.134E-I-05 


2.30 


100.28 


1122 


c 


1.20E+02 


0.83628 


6.9999E-09 


21.1155 


2.22E-06 


1.115E+05 


1.75 


98.63 


1123 


c 


1.12E+02 


0.82528 


6.7743E-09 


21.1352 


2.09E-06 


1.096E+05 


1.32 


96.99 


1124 


c 


1.04E+02 


0.81748 


6.5834E-09 


21.1566 


1.96E-06 


1.080E+05 


1.00 


95.77 


1125 


c 


l.OOE+02 


0.81258 


6.4480E-09 


21.1676 


1.89E-06 


1.069E+05 


0.80 


94.97 


1126 


c 


9.62E+01 


0.80607 


6.2478E-09 


21.1788 


1.82E-06 


1.054E-I-05 


0.67 


93.86 


1127 


c 


8.92E+01 


0.80011 


6.0434E-09 


21.1996 


1.69E-06 


1.038E+05 


0.52 


92.80 


1128 


c 


8.57E+01 


0.79626 


5.9007E-09 


21.2104 


1.63E-06 


1.028E+05 


0.42 


92.09 


1129 


c 


8.22E+01 


0.79196 


5.7334E-09 


21.2214 


1.57E-06 


1.016E+05 


0.35 


91.30 


1130 


c 


7.79E+01 


0.78715 


5.5390E-09 


21.2356 


1.49E-06 


1.002E+05 


0.27 


90.43 


1131 


c 


7.36E+01 


0.78269 


5.3475E-09 


21.2503 


1.41E-06 


9.880E-I-04 


0.21 


89.60 


1132 


c 


6.98E+01 


0.77884 


5.1603E-09 


21.2638 


1.34E-06 


9.747E+04 


0.15 


88.83 


1133 


c 


6.61E+01 


0.77531 


4.9665E-09 


21.2778 


1.27E-06 


9.609E+04 


0.11 


88.05 


1134 


c 


6.23E+01 


0.77266 


4.8046E-09 


21.2923 


1.20E-06 


9.494E+04 


0.07 


87.42 


1135 


c 


6.00E-1-01 


0.77075 


4.6775E-09 


21.3017 


1.16E-06 


9.403E-I-04 


0.05 


86.93 


1136 


c 


5.77E+01 


0.76804 


4.4821E-09 


21.3112 


1.12E-06 


9.261E-I-04 


0.03 


86.19 


1137 


c 


5.32E+01 


0.76529 


4.2679E-09 


21.3311 


1.03E-06 


9.105E+04 


0.01 


85.39 


1138 


R 


5.04E+01 


0.76325 


4.1031E-09 


21.3436 


9.78E-07 


8.983E+04 


0.00 


84.78 


1139 


R 


4.77E-I-01 


0.76129 


3.9433E-09 


21.3566 


9.26E-07 


8.863E+04 


0.00 


84.20 


1140 


R 


4.53E+01 


0.75959 


3.7889E-09 


21.3687 


8.80E-07 


8.746E+04 


0.00 


83.64 


1141 


R 


4.29E+01 


0.75813 


3.6339E-09 


21.3813 


8.33E-07 


8.627E+04 


0.00 


83.08 


1142 


R 


4.06E+01 


0.75688 


3.4786E-09 


21.3939 


7.89E-07 


8.505E+04 


0.00 


82.52 


1143 


R 


3.83E+01 


0.75581 


3.3202E-09 


21.4071 


7.44E-07 


8.378E+04 


0.00 


81.95 


1144 


R 


3.60E+01 


0.75487 


3.1587E-09 


21.4209 


6.99E-07 


8.246E-I-04 


0.00 


81.37 


1145 


R 


3.37E+01 


0.75414 


3.0163E-09 


21.4353 


6.54E-07 


8.126E+04 


0.00 


80.85 


1146 


R 


3.20E+01 


0.75353 


2.8941E-09 


21.4464 


6.22E-07 


8.021E+04 


0.00 


80.40 


1147 


R 


3.04E+01 


0.75290 


2.7700E-09 


21.4579 


5.89E-07 


7.913E+04 


0.00 


79.94 


1148 


R 


2.87E+01 


0.75225 


2.6285E-09 


21.4699 


5.56E-07 


7.786E+04 


0.00 


79.41 


1149 


R 


2.67E+01 


0.75167 


2.4687E-09 


21.4856 


5.16E-07 


7.638E-I-04 


0.00 


78.80 


1150 


R 


2.46E+01 


0.75132 


2.3386E-09 


21.5023 


4.75E-07 


7.514E+04 


0.00 


78.31 


1151 


R 


2.34E+01 


0.75114 


2.2398E-09 


21.5128 


4.51E-07 


7.417E+04 


0.00 


77.93 


1152 


R 


2.22E+01 


0.75102 


2.1310E-09 


21.5238 


4.27E-07 


7.308E+04 


0.00 


77.51 


1153 


R 


2.08E+01 


0.75100 


2.0117E-09 


21.5373 


3.98E-07 


7.185E+04 


0.00 


77.04 


1154 


R 


1.93E+01 


0.75107 


1.9108E-09 


21.5515 


3.70E-07 


7.077E-I-04 


0.00 


76.64 


1155 


R 


1.84E+01 


0.75118 


1.8301E-09 


21.5615 


3.51E-07 


6.989E+04 


0.00 


76.32 


1156 


R 


1.75E+01 


0.75135 


1.7497E-09 


21.5717 


3.33E-07 


6.899E+04 


0.00 


75.99 


1157 


R 


1.65E+01 


0.75159 


1.6682E-09 


21.5823 


3.14E-07 


6.805E+04 


0.00 


75.66 


1158 


R 


1.56E-1-01 


0.75193 


1.5858E-09 


21.5934 


2.96E-07 


6.707E-I-04 


0.00 


75.32 


1159 


R 


1.47E+01 


0.75232 


1.5107E-09 


21.6051 


2.77E-07 


6.616E-I-04 


0.00 


75.02 


1160 


R 


1.39E+01 


0.75277 


1.4434E-09 


21.6148 


2.63E-07 


6.531E+04 


0.00 


74.74 


1161 


R 


1.32E+01 


0.75333 


1.3753E-09 


21.6250 


2.48E-07 


6.444E+04 


0.00 


74.45 


1162 


R 


1.25E+01 


0.75402 


1.3065E-09 


21.6358 


2.33E-07 


6.352E+04 


0.00 


74.16 


1163 


R 


1.17E+01 


0.75502 


1.2290E-09 


21.6470 


2.19E-07 


6.246E+04 


0.00 


73.83 


1164 


R 


1.09E+01 


0.75652 


1.1428E-09 


21.6617 


2.01E-07 


6.122E+04 


0.00 


73.45 


1165 


R 


9.97E+00 


0.75795 


1.0761E-09 


21.6774 


1.83E-07 


6.023E+04 


0.00 


73.15 


1166 


R 


9.49E+00 


0.75910 


1.0296E-09 


21.6863 


1.73E-07 


5.950E+04 


0.00 


72.94 


1167 


R 


9.01E+00 


0.76037 


9.8282E-10 


21.6955 


1.64E-07 


5.876E+04 


0.00 


72.72 


1168 


R 


8.54E+00 


0.76193 


9.3051E-10 


21.7052 


1.54E-07 


5.789E+04 


0.00 


72.46 


1169 


R 


7.96E+00 


0.76379 


8.7258E-10 


21.7176 


1.42E-07 


5.689E+04 


0.00 


72.17 


1170 


R 


7.38E+00 


0.76562 


8.1920E-10 


21.7309 


1.31E-07 


5.593E+04 


0.00 


71.88 


1171 


R 


6.89E+00 


0.76732 


7.7565E-10 


21.7426 


1.21E-07 


5.511E+04 


0.00 


71.64 


1172 


R 


6.51E+00 


0.76905 


7.3709E-10 


21.7523 


1.14E-07 


5.435E-I-04 


0.00 


71.41 


1173 


R 


6.12E+00 


0.77070 


7.0359E-10 


21.7626 


1.06E-07 


5.367E+04 


0.00 


71.21 


1174 


R 


5.84E+00 


0.77221 


6.7523E-10 


21.7704 


l.OOE-07 


5.308E+04 


0.00 


71.03 
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Grid point 


STAT 


T 


X 


P 






M. - Mr [g] 


T[K] 


Vc [km s"^] 


Cs [km s~^] 


1175 


R 


5.56E+00 


0.77414 


6.4136E- 


10 


21.7785 


9.48E-08 


5.234E+04 


0.00 


70.80 


1176 


R 


5.17E+00 


0.77624 


6.0673E- 


10 


21.7904 


8.71E-08 


5.156E+04 


0.00 


70.55 


1177 


R 


4.87E-I-00 


0.77813 


5.7676E- 


10 


21.7999 


8.13E-08 


5.085E+04 


0.00 


70.32 


1178 


R 


4.58E+00 


0.77975 


5.5149E- 


10 


21.8098 


7.54E-08 


5.024E+04 


0.00 


70.11 


1179 


R 


4.37E+00 


0.78112 


5.3095E- 


10 


21.8169 


7.15E-08 


4.972E+04 


0.00 


69.94 


1180 


R 


4.17E+00 


0.78266 


5.1037E- 


10 


21.8243 


6.75E-08 


4.918E+04 


0.00 


69.75 


1181 


R 


3.96E+00 


0.78445 


4.8907E- 


10 


21.8320 


6.35E-08 


4.861E+04 


0.00 


69.55 


1182 


R 


3.74E+00 


0.78647 


4.6708E- 


10 


21.8405 


5.93E-08 


4.799E-I-04 


0.00 


69.32 


1183 


R 


3.53E+00 


0.78864 


4.4509E- 


10 


21.8494 


5.51E-08 


4.735E+04 


0.00 


69.07 


1184 


R 


3.31E+00 


0.79092 


4.231 1E- 


10 


21.8587 


5.08E-08 


4.669E+04 


0.00 


68.80 


1185 


R 


3.09E+00 


0.79285 


4.0457E- 


10 


21.8686 


4.66E-08 


4.610E+04 


0.00 


68.55 


1186 


R 


2.94E+00 


0.79437 


3.8946E- 


10 


21.8756 


4.37E-08 


4.560E+04 


0.00 


68.33 


1187 


R 


2.78E+00 


0.79580 


3.7432E- 


10 


21.8830 


4.08E-08 


4.509E-I-04 


0.00 


68.10 


1188 


R 


2.63E+00 


0.79704 


3.5915E- 


10 


21.8906 


3.79E-08 


4.456E+04 


0.00 


67.85 


1189 


R 


2.48E+00 


0.79808 


3.4600E- 


10 


21.8985 


3.49E-08 


4.409E+04 


0.00 


67.62 


1190 


R 


2.37E+00 


0.79898 


3.3489E- 


10 


21.9045 


3.28E-08 


4.367E+04 


0.00 


67.42 


1191 


R 


2.26E+00 


0.80002 


3.2165E- 


10 


21.9107 


3.07E-08 


4.317E+04 


0.00 


67.16 


1192 


R 


2.10E+00 


0.80087 


3.0983E- 


10 


21.9195 


2.78E-08 


4.271E-I-04 


0.00 


66.93 


1193 


R 


2.02E+00 


0.80139 


3.0157E- 


10 


21.9244 


2.62E-08 


4.238E+04 


0.00 


66.76 


1194 


R 


1.94E+00 


0.80182 


2.9326E- 


10 


21.9295 


2.47E-08 


4.204E+04 


0.00 


66.58 


1195 


R 


1.86E+00 


0.80213 


2.849 1E- 


10 


21.9346 


2.31E-08 


4.169E+04 


0.00 


66.40 


1196 


R 


1.78E-1-00 


0.80227 


2.7606E- 


10 


21.9400 


2.16E-08 


4.131E-I-04 


0.00 


66.20 


1197 


R 


1.69E+00 


0.80215 


2.6667E- 


10 


21.9460 


1.99E-08 


4.091E-I-04 


0.00 


65.99 


1198 


R 


1.59E+00 


0.80180 


2.5928E- 


10 


21.9523 


1.81E-08 


4.058E+04 


0.00 


65.83 


1199 


R 


1.54E+00 


0.80112 


2.5125E- 


10 


21.9559 


1.72E-08 


4.023E+04 


0.00 


65.65 


1200 


R 


1.44E+00 


0.80003 


2.4310E- 


10 


21.9634 


1.53E-08 


3.986E+04 


0.00 


65.48 


1201 


R 


1.39E+00 


0.79900 


2.3729E- 


10 


21.9672 


1.43E-08 


3.960E+04 


0.00 


65.36 


1202 


R 


1.34E+00 


0.79775 


2.3130E- 


10 


21.9715 


1.32E-08 


3.933E+04 


0.00 


65.25 


1203 


R 


1.29E+00 


0.79556 


2.2262E- 


10 


21.9757 


1.22E-08 


3.894E+04 


0.00 


65.10 


1204 


R 


1.19E+00 


0.79285 


2.1391E- 


10 


21.9841 


1.03E-08 


3.854E+04 


0.00 


64.98 


1205 


R 


1.14E+00 


0.78970 


2.0556E- 


10 


21.9883 


9.40E-09 


3.817E-I-04 


0.00 


64.90 


1206 


R 


1.05E+00 


0.78603 


1.9736E- 


10 


21.9964 


7.68E-09 


3.780E+04 


0.00 


64.86 


1207 


R 


l.OlE+00 


0.78257 


1.9068E- 


10 


22.0003 


6.88E-09 


3.751E+04 


0.00 


64.86 


1208 


R 


9.50E-01 


0.77882 


1.8424E- 


10 


22.0064 


5.69E-09 


3.723E+04 


0.00 


64.89 


1209 


R 


9.17E-01 


0.77546 


1.7903E- 


10 


22.0098 


5.03E-09 


3.700E+04 


0.00 


64.95 


1210 


R 


8.76E-01 


0.77201 


1.7411E- 


10 


22.0141 


4.24E-09 


3.680E-I-04 


0.00 


65.03 


1211 


R 


8.49E-01 


0.76795 


1.6877E- 


10 


22.0172 


3.69E-09 


3.658E+04 


0.00 


65.15 


1212 


R 


8.06E-01 


0.76388 


1.6380E- 


10 


22.0221 


2.83E-09 


3.637E+04 


0.00 


65.30 


1213 


R 


7.85E-01 


76046 


1.5988E- 


10 


22.0245 


2.42E-09 


3.622E+04 


0.00 


65.44 


1214 


R 


7.57E-01 


0.75691 


1.5601E- 


10 


22.0280 


1.85E-09 


3.606E+04 


0.00 


65.60 


1215 


R 


7.39E-01 


0.75246 


1.5144E- 


10 


22.0303 


1.47E-09 


3.589E-I-04 


0.00 


65.83 


1216 


R 


7.04E-01 


0.74810 


1.4719E- 


10 


22.0348 


7.64E-10 


3.573E+04 


0.00 


66.07 


1217 


R 


6.91E-01 


0.74443 


1.4378E- 


10 


22.0366 


4.96E-10 


3.560E+04 


0.00 


66.29 


1218 


R 


6.67E-01 


0.74192 


1.4152E- 


10 


22.0399 


O.OOE+00 


3.552E+04 


0.00 


66.45 
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Table A.2. Outermost 1 Rq of a 20 Mq model at solar metallicity (Z=0.02). 
The table shows the physical variables at 1=6.41x10'' years. Columns contain the 
progressive grid point of the model, the status (R=radiative, C=convective), the 
optical depth t, the opacity x, the density p, the radius R, the value M, - 
(where M, is the total stellar mass and Mr is the mass coordinate), the tempera- 
ture T, the convective velocity and the local sound speed Cj. All the values are 
in cgs units if not otherwise specified. 



id point 


o lAl 




T 








P 




K [RqJ 


A4 


A/I frrl 

- Mr LgJ 




1[KJ 


V(. [km s ] 


Cj [km s 


1009 


R 


1, 


.88E+04 


1 


.15006 


9 


.1148E- 


■07 


9.4516 


3, 


.91E-05 


3, 


.415E+05 


0.00 


110.14 


1010 


R 


1, 


.80E+04 


1, 


.16843 


8, 


.467 IE- 


■07 


9.4634 


3, 


.70E-05 


3, 


.363E+05 


0.00 


109.93 


1011 


R 


1 


.68E-H04 


1 


.18658 


7, 


.8819E- 


■07 


9.4826 


3, 


.39E-05 


3, 


.313E-I-05 


0.00 


109.74 


1012 


R 


1 


.61E-I-04 


1 


.19963 


7, 


.4912E- 


■07 


9.4929 


3, 


.23E-05 


3, 


.279E-I-05 


0.00 


109.61 


1013 


R 


1, 


.55E+04 


1, 


.21788 


6 


.9847E- 


■07 


9.5037 


3, 


.08E-05 


3, 


.232E+05 


0.00 


109.45 


1014 


R 


1, 


.45E+04 


1, 


.24209 


6 


.3760E- 


■07 


9.5220 


2, 


.83E-05 


3, 


.172E+05 


0.00 


109.26 


1015 


R 


1, 


.35E+04 


1, 


.26170 


5, 


.9299E- 


■07 


9.5413 


2, 


.59E-05 


3, 


.125E+05 


0.00 


109.14 


1016 


R 


1 


.30E-H04 


1 


.27584 


5 


.6324E- 


■07 


9.5516 


2, 


.47E-05 


3, 


.093E-I-05 


0.00 


109.06 


1017 


R 


1 


.24E-I-04 


1 


.29061 


5 


.3416E- 


■07 


9.5624 


2, 


.36E-05 


3, 


.060E-I-05 


0.00 


108.98 


1018 


R 


1 


.19E-H04 


1, 


.30554 


5, 


.0665E- 


■07 


9.5733 


2, 


.24E-05 


3, 


.027E+05 


0.00 


108.92 


1019 


R 


1, 


.15E+04 


1, 


.32091 


4, 


.8003E- 


■07 


9.5841 


2, 


.14E-05 


2, 


.994E+05 


0.00 


108.86 


1020 


R 


1, 


.lOE+04 


1, 


.33683 


4, 


.5439E- 


■07 


9.5954 


2, 


.03E-05 


2, 


.961E+05 


0.00 


108.81 


1021 


R 


1 


.05E-H04 


1 


.35328 


4 


.2968E- 


■07 


9.6065 


1, 


.93E-05 


2, 


.928EH-05 


0.00 


108.77 


1022 


R 


1 


.OlE-i-04 


1 


.37094 


4 


.0503E- 


■07 


9.6182 


1, 


.83E-05 


2, 


.893E-I-05 


0.00 


108.73 


1023 


R 


9 


.62E-H03 


1, 


.38833 


3, 


.8248E- 


■07 


9.6306 


1, 


.73E-05 


2, 


.860E-I-05 


0.00 


108.71 


1024 


R 


9 


.24E-H03 


1, 


.40532 


3, 


.6200E- 


■07 


9.6416 


1, 


.65E-05 


2, 


.828E-H05 


0.00 


108.69 


1025 


R 


8, 


.85E-H03 


1, 


.42367 


3, 


.4159E- 


■07 


9.6531 


1, 


.57E-05 


2, 


.795E-H05 


0.00 


108.68 


1026 


R 


8 


.45E-I-03 


1 


.44389 


3 


.2125E- 


■07 


9.6653 


1, 


.48E-05 


2, 


.761EH-05 


0.00 


108.67 


1027 


R 


8 


.05E-I-03 


1 


.46417 


3 


.0276E- 


■07 


9.6783 


1, 


.40E-05 


2, 


.728E-I-05 


0.00 


108.67 


1028 


R 


7, 


.72E+03 


1, 


.48419 


2, 


.861 IE- 


■07 


9.6896 


1, 


.33E-05 


2, 


.697E+05 


0.00 


108.68 


1029 


R 


7, 


.38E-I-03 


1 


.50595 


2, 


.6954E- 


■07 


9.7016 


1, 


.26E-05 


2, 


.664E-I-05 


0.00 


108.68 


1030 


R 


7, 


.04E-I-03 


1 


.52815 


2, 


.5408E- 


■07 


9.7142 


1, 


.19E-05 


2, 


.632E-I-05 


0.00 


108.68 


1031 


R 


6 


.73E-H03 


1 


.55060 


2, 


.3969E- 


■07 


9.7260 


1, 


.13E-05 


2, 


.601E-I-05 


0.00 


108.68 


1032 


R 


6 


.42E+03 


1, 


.57310 


2, 


.2639E- 


■07 


9.7385 


1, 


.07E-05 


2, 


.570E+05 


0.00 


108.67 


1033 


R 


6 


.15E+03 


1, 


.59542 


2, 


.1415E- 


■07 


9.7499 


1, 


.02E-05 


2, 


.541E+05 


0.00 


108.66 


1034 


R 


5, 


.88E-H03 


1 


.62281 


2, 


.0030E- 


■07 


9.7619 


9, 


.69E-06 


2, 


.505E-I-05 


0.00 


108.63 


1035 


R 


5, 


.53E-I-03 


1 


.65803 


1 


.8493E- 


■07 


9.7784 


9, 


.02E-06 


2, 


.463E-I-05 


0.00 


108.57 


1036 


R 


5 


.17E+03 


1 


.69089 


1, 


.7276E- 


■07 


9.7961 


8, 


.35E-06 


2, 


.426E+05 


0.00 


108.48 


1037 


R 


4, 


.94E+03 


1, 


.72391 


1, 


.621 IE- 


■07 


9.8077 


7, 


.94E-06 


2, 


.392E+05 


0.00 


108.36 


1038 


R 


4, 


.64E+03 


1, 


.76662 


1, 


.5016E- 


■07 


9.8242 


7, 


.39E-06 


2, 


.350E+05 


0.00 


108.14 


1039 


R 


4, 


.32E-H03 


1 


.80536 


1 


.4067E- 


■07 


9.8420 


6, 


.84E-06 


2, 


.314E-I-05 


0.00 


107.88 


1040 


R 


4 


.12E-I-03 


1 


.83771 


1 


.3352E- 


■07 


9.8539 


6, 


.50E-06 


2, 


.284E-I-05 


0.00 


107.60 


1041 


R 


3, 


.91E+03 


1, 


.86954 


1, 


.2702E- 


■07 


9.8664 


6, 


.15E-06 


2, 


.255E+05 


0.00 


107.28 


1042 


R 


3, 


.73E+03 


1, 


.90044 


1, 


.2115E- 


■07 


9.8776 


5, 


.85E-06 


2, 


.226E+05 


0.00 


106.90 


1043 


R 


3, 


.55E-H03 


1, 


.93471 


1, 


.1542E- 


■07 


9.8894 


5, 


.56E-06 


2, 


.197E-f05 


0.00 


106.44 


1044 


R 


3 


.37E-H03 


1 


.97305 


1 


.0986E- 


■07 


9.9017 


5, 


.26E-06 


2, 


.165E-I-05 


0.00 


105.86 


1045 


R 


3, 


.18E-I-03 


2, 


.01222 


1 


.0492E- 


■07 


9.9146 


4, 


.96E-06 


2, 


.134E-I-05 


0.00 


105.21 


1046 


R 


3, 


.01E-H03 


2, 


.05131 


1, 


.0056E- 


■07 


9.9260 


4, 


.71E-06 


2, 


.105E+05 


0.00 


104.47 


1047 


C 


2, 


.85E-I-03 


2, 


.09321 


9 


.6390E- 


■08 


9.9379 


4, 


.46E-06 


2, 


.073EH-05 


0.27 


103.59 


1048 


C 


2, 


.68E-H03 


2, 


.13709 


9 


.2418E- 


■08 


9.9502 


4, 


.21E-06 


2, 


.039E-I-05 


3.48 


102.53 


1049 


C 


2, 


.50E-I-03 


2, 


.17706 


8 


.8976E- 


■08 


9.9630 


3. 


.96E-06 


2. 


.007E-I-05 


5.74 


101.40 


1050 


C 


2, 


.36E-I-03 


2, 


.21180 


8 


.6007E- 


■08 


9.9739 


3, 


.76E-06 


1, 


.976E-I-05 


7.30 


100.24 


1051 


C 


2, 


.21E+03 


2, 


.24604 


8, 


.3153E- 


■08 


9.9851 


3, 


.55E-06 


1, 


.943E+05 


8.51 


98.93 


1052 


C 


2, 


.07E-H03 


2, 


.27681 


8 


.0595E- 


■08 


9.9967 


3, 


.35E-06 


1, 


.911E-I-05 


9.43 


97.55 


1053 


C 


1 


.94E-I-03 


2, 


.30324 


7, 


.8309E- 


■08 


10.0070 


3, 


.17E-06 


1, 


.879E-I-05 


10.10 


96.14 


1054 


c 


1, 


.80E+03 


2, 


.32459 


7, 


.6279E- 


■08 


10.0177 


2, 


.99E-06 


1, 


.848E+05 


10.53 


94.69 


1055 


c 


1, 


.69E+03 


2, 


.34028 


7, 


.4483E- 


■08 


10.0269 


2, 


.84E-06 


1, 


.818E+05 


10.74 


93.24 


1056 


c 


1, 


.58E+03 


2, 


.35064 


7, 


.2749E- 


■08 


10.0363 


2, 


.69E-06 


1, 


.786E+05 


10.81 


91.67 


1057 


c 


1 


.47E-H03 


2, 


.35334 


7, 


.1321E- 


■08 


10.0460 


2, 


.54E-06 


1, 


.758E-I-05 


10.61 


90.24 


1058 


c 


1 


.39E-I-03 


2, 


.34958 


7, 


.0171E- 


■08 


10.0528 


2, 


.43E-06 


1, 


.734E-I-05 


10.26 


89.00 


1059 


c 


1 


.31E-H03 


2, 


.33995 


6 


.9027E- 


■08 


10.0598 


2, 


.33E-06 


1, 


.708E+05 


9.87 


87.70 


1060 


c 


1 


.23E-H03 


2, 


.32417 


6 


.7878E- 


■08 


10.0668 


2, 


.22E-06 


1, 


.682E+05 


9.34 


86.34 


1061 


c 


1, 


.16E-H03 


2, 


.30414 


6 


.6802E- 


■08 


10.0740 


2, 


.12E-06 


1, 


.657E+05 


8.69 


85.05 


1062 


c 


1 


.09E-H03 


2, 


.28134 


6 


.579 lE- 


■08 


10.0801 


2, 


.03E-06 


1, 


.633E-I-05 


8.00 


83.83 
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Grid point 


STAT 


T 


X 


P 




- M, [g] 


T[K] 


Vc [km s"^] 


Cs [km s ^] 


1063 


c 


1.03E+03 


2.25516 


6.4766E-08 


10.0863 


1.94E-06 


1.610E+05 


7.32 


82.61 


1064 


c 


9.71E+02 


2.22633 


6.3719E-08 


10.0923 


1.86E-06 


1.585E+05 


6.62 


81.39 


1065 


c 


9.12E-I-02 


2.19204 


6.2614E-08 


10.0984 


1.77E-06 


1.561E+05 


5.91 


80.16 


1066 


c 


8.55E+02 


2.15002 


6.1433E-08 


10.1046 


1.69E-06 


1.535E+05 


5.16 


78.90 


1067 


c 


8.00E+02 


2.10271 


6.0183E-08 


10.1109 


1.60E-06 


1.509E+05 


4.39 


77.68 


1068 


c 


7.48E+02 


2.05208 


5.8853E-08 


10.1170 


1.52E-06 


1.483E+05 


3.66 


76.48 


1069 


c 


6.97E+02 


2.01013 


5.7714E-08 


10.1233 


1.44E-06 


1.462E+05 


3.03 


75.53 


1070 


c 


6.68E+02 


1.97885 


5.6821E-08 


10.1270 


1.40E-06 


1.447E-I-05 


2.59 


74.84 


1071 


c 


6.40E+02 


1.94778 


5.5882E-08 


10.1308 


1.35E-06 


1.431E+05 


2.26 


74.14 


1072 


c 


6.11E+02 


1.91749 


5.4897E-08 


10.1346 


1.30E-06 


1.415E+05 


1.97 


73.45 


1073 


c 


5.84E+02 


1.88866 


5.3907E-08 


10.1386 


1.25E-06 


1.399E+05 


1.71 


72.79 


1074 


c 


5.59E+02 


1.86001 


5.2917E-08 


10.1422 


1.21E-06 


1.383E+05 


1.49 


72.15 


1075 


c 


5.34E+02 


1.83056 


5.1882E-08 


10.1459 


1.17E-06 


1.368E-I-05 


1.28 


71.52 


1076 


c 


5.10E+02 


1.80057 


5.0800E-08 


10.1497 


1.13E-06 


1.352E+05 


1.09 


70.88 


1077 


c 


4.86E+02 


1.76619 


4.9509E-08 


10.1536 


1.08E-06 


1.333E+05 


0.91 


70.17 


1078 


c 


4.57E+02 


1.72789 


4.7985E-08 


10.1586 


1.03E-06 


1.312E+05 


0.73 


69.36 


1079 


c 


4.29E+02 


1.69264 


4.6470E-08 


10.1638 


9.72E-07 


1.292E+05 


0.56 


68.61 


1080 


c 


4.04E+02 


1.66094 


4.4977E-08 


10.1686 


9.24E-07 


1.272E-I-05 


0.43 


67.90 


1081 


c 


3.80E+02 


1.63122 


4.3415E-08 


10.1735 


8.75E-07 


1.252E+05 


0.32 


67.18 


1082 


c 


3.56E+02 


1.60288 


4.1784E-08 


10.1786 


8.26E-07 


1.232E+05 


0.23 


66.46 


1083 


c 


3.33E+02 


1.57488 


4.0035E-08 


10.1838 


7.78E-07 


1.210E+05 


0.16 


65.72 


1084 


c 


3.09E-1-02 


1.54737 


3.8160E-08 


10.1896 


7.26E-07 


1.186E-I-05 


0.10 


64.94 


1085 


c 


2.86E+02 


1.52243 


3.6302E-08 


10.1957 


6.75E-07 


1.163E-I-05 


0.05 


64.18 


1086 


c 


2.65E+02 


1.49992 


3.4470E-08 


10.2015 


6.29E-07 


1.140E+05 


0.01 


63.44 


1087 


R 


2.45E+02 


1.48240 


3.2925E-08 


10.2075 


5.83E-07 


1.121E+05 


0.00 


62.83 


1088 


R 


2.32E+02 


1.46981 


3.1697E-08 


10.2115 


5.54E-07 


1.105E+05 


0.00 


62.35 


1089 


R 


2.19E+02 


1.45895 


3.0510E-08 


10.2156 


5.24E-07 


1.090E+05 


0.00 


61.88 


1090 


R 


2.08E+02 


1.44955 


2.9369E-08 


10.2194 


4.99E-07 


1.075E+05 


0.00 


61.43 


1091 


R 


1.97E+02 


1.44085 


2.8205E-08 


10.2234 


4.73E-07 


1.060E+05 


0.00 


60.97 


1092 


R 


1.86E+02 


1.43173 


2.6869E-08 


10.2275 


4.47E-07 


1.042E+05 


0.00 


60.43 


1093 


R 


1.72E+02 


1.42204 


2.5352E-08 


10.2329 


4.15E-07 


1.022E+05 


0.00 


59.81 


1094 


R 


1.59E+02 


1.41419 


2.4109E-08 


10.2385 


3.82E-07 


1.005E+05 


0.00 


59.29 


1095 


R 


1.50E+02 


1.40808 


2.3160E-08 


10.2421 


3.63E-07 


9.911E+04 


0.00 


58.89 


1096 


R 


1.42E+02 


1.40244 


2.2244E-08 


10.2458 


3.44E-07 


9.779E+04 


0.00 


58.50 


1097 


R 


1.35E+02 


1.39730 


2.1365E-08 


10.2493 


3.27E-07 


9.650E+04 


0.00 


58.12 


1098 


R 


1.28E+02 


1.39246 


2.0485E-08 


10.2529 


3.10E-07 


9.519E-I-04 


0.00 


57.73 


1099 


R 


1.21E+02 


1.38794 


1.9605E-08 


10.2565 


2.93E-07 


9.385E+04 


0.00 


57.34 


1100 


R 


1.14E+02 


1.38372 


1.8710E-08 


10.2603 


2.76E-07 


9.246E+04 


0.00 


56.93 


1101 


R 


1.08E+02 


1.37985 


1.7798E-08 


10.2642 


2.60E-07 


9.101E+04 


0.00 


56.51 


1102 


R 


l.OlE+02 


1.37684 


1.6995E-08 


10.2684 


2.43E-07 


8.970E+04 


0.00 


56.13 


1103 


R 


9.59E+01 


1.37462 


1.6306E-08 


10.2716 


2.31E-07 


8.855E-I-04 


0.00 


55.80 


1104 


R 


9.10E+01 


1.37291 


1.5607E-08 


10.2749 


2.19E-07 


8.736E+04 


0.00 


55.46 


1105 


R 


8.61E+01 


1.37164 


1.481 lE-08 


10.2783 


2.07E-07 


8.597E+04 


0.00 


55.07 


1106 


R 


8.01E+01 


1.37106 


1.3914E-08 


10.2829 


1.92E-07 


8.435E+04 


0.00 


54.60 


1107 


R 


7.41E-I-01 


1.37119 


1.3184E-08 


10.2877 


1.77E-07 


8.299E-I-04 


0.00 


54.22 


1108 


R 


7.04E+01 


1.37161 


1.2631E-08 


10.2907 


1.68E-07 


8.192E-I-04 


0.00 


53.91 


1109 


R 


6.68E+01 


1.37234 


1.2022E-08 


10.2938 


1.59E-07 


8.072E+04 


0.00 


53.57 


1110 


R 


6.26E+01 


1.37335 


1.1354E-08 


10.2977 


1.48E-07 


7.936E+04 


0.00 


53.18 


nil 


R 


5.83E+01 


1.37444 


1.0789E-08 


10.3018 


1.37E-07 


7.817E+04 


0.00 


52.85 


1112 


R 


5.55E+01 


1.37559 


1.0338E-08 


10.3047 


1.30E-07 


7.720E+04 


0.00 


52.57 


1113 


R 


5.27E+01 


1.37700 


9.8872E-09 


10.3076 


1.24E-07 


7.620E+04 


0.00 


52.29 


1114 


R 


5.00E+01 


1.37874 


9.4310E-09 


10.3107 


1.17E-07 


7.516E+04 


0.00 


51.99 


1115 


R 


4.72E+01 


1.38088 


8.9689E-09 


10.3139 


l.lOE-07 


7.407E+04 


0.00 


51.69 


1116 


R 


4.44E+01 


1.38319 


8.5480E-09 


10.3173 


1.03E-07 


7.305E-I-04 


0.00 


51.40 


1117 


R 


4.22E+01 


1.38562 


8.1697E-09 


10.3201 


9.76E-08 


7.211E+04 


0.00 


51.13 


1118 


R 


4.00E+01 


1.38845 


7.7869E-09 


10.3230 


9.22E-08 


7.112E+04 


0.00 


50.85 


1119 


R 


3.78E+01 


1.39177 


7.3995E-09 


10.3261 


8.67E-08 


7.010E+04 


0.00 


50.56 


1120 


R 


3.56E+01 


1.39617 


6.9626E-09 


10.3293 


8.12E-08 


6.890E+04 


0.00 


50.23 


1121 


R 


3.29E+01 


1.40204 


6.4740E-09 


10.3336 


7.45E-08 


6.750E-I-04 


0.00 


49.83 


1122 


R 


3.02E+01 


1.40739 


6.0948E-09 


10.3381 


6.79E-08 


6.636E+04 


0.00 


49.51 


1123 


R 


2.88E+01 


1.41162 


5.8291E-09 


10.3406 


6.43E-08 


6.554E+04 


0.00 


49.27 
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Grid point 


STAT 


T 


X 


P 




M. - M, [g] 


T[K] 


Vc [km s"^] 


Cs [km s ^] 


1124 


R 


2.73E+01 


1.41634 


5.5609E-09 


10.3433 


6.08E-08 


6.468E+04 


0.00 


49.03 


1125 


R 


2.59E+01 


1.42226 


5.2595E-09 


10.3461 


5.72E-08 


6.368E+04 


0.00 


48.74 


1126 


R 


2.41E+01 


1.42985 


4.9241E-09 


10.3497 


5.29E-08 


6.252E+04 


0.00 


48.41 


1127 


R 


2.23E+01 


1.43837 


4.6133E-09 


10.3536 


4.86E-08 


6.140E+04 


0.00 


48.09 


1128 


R 


2.08E+01 


1.44656 


4.3587E-09 


10.3570 


4.50E-08 


6.043E+04 


0.00 


47.81 


1129 


R 


1.96E+01 


1.45480 


4.1322E-09 


10.3598 


4.22E-08 


5.954E+04 


0.00 


47.54 


1130 


R 


1.84E+01 


1.46274 


3.9346E-09 


10.3628 


3.93E-08 


5.873E+04 


0.00 


47.30 


1131 


R 


1.76E+01 


1.47005 


3.7666E-09 


10.3651 


3.73E-08 


5.802E-I-04 


0.00 


47.08 


1132 


R 


1.67E+01 


1.47948 


3.5652E-09 


10.3675 


3.52E-08 


5.712E+04 


0.00 


46.80 


1133 


R 


1.55E+01 


1.48989 


3.3581E-09 


10.3710 


3.24E-08 


5.617E+04 


0.00 


46.50 


1134 


R 


1.45E+01 


1.50018 


3.1779E-09 


10.3738 


3.02E-08 


5.529E+04 


0.00 


46.22 


1135 


R 


1.36E+01 


1.51024 


3.0255E-09 


10.3767 


2.80E-08 


5.452E+04 


0.00 


45.96 


1136 


R 


1.29E+01 


1.51938 


2.9012E-09 


10.3789 


2.66E-08 


5.386E-I-04 


0.00 


45.74 


1137 


R 


1.23E+01 


1.52942 


2.7765E-09 


10.3811 


2.51E-08 


5.318E+04 


0.00 


45.51 


1138 


R 


1.16E+01 


1.54073 


2.6471E-09 


10.3834 


2.36E-08 


5.244E+04 


0.00 


45.24 


1139 


R 


1.09E+01 


1.55335 


2.5130E-09 


10.3860 


2.20E-08 


5.164E+04 


0.00 


44.94 


1140 


R 


1.02E+01 


1.56677 


2.3784E-09 


10.3887 


2.05E-08 


5.079E+04 


0.00 


44.61 


1141 


R 


9.50E+00 


1.58081 


2.2433E-09 


10.3916 


1.89E-08 


4.988E-I-04 


0.00 


44.25 


1142 


R 


8.78E+00 


1.59453 


2.1289E-09 


10.3946 


1.73E-08 


4.908E+04 


0.00 


43.91 


1143 


R 


8.28E+00 


1.60744 


2.0356E-09 


10.3968 


1.62E-08 


4.838E+04 


0.00 


43.60 


1144 


R 


7.77E+00 


1.62179 


1.9422E-09 


10.3991 


1.51E-08 


4.765E+04 


0.00 


43.25 


1145 


R 


7.26E-1-00 


1.63733 


1.8489E-09 


10.4016 


1.41E-08 


4.688E-I-04 


0.00 


42.87 


1146 


R 


6.75E+00 


1.65130 


1.7685E-09 


10.4041 


1.30E-08 


4.618E-I-04 


0.00 


42.50 


1147 


R 


6.37E+00 


1.66300 


1.7008E-09 


10.4060 


1.22E-08 


4.555E+04 


0.00 


42.15 


1148 


R 


5.99E+00 


1.67593 


1.6205E-09 


10.4080 


1.14E-08 


4.477E+04 


0.00 


41.69 


1149 


c 


5.46E-I-00 


1.68633 


1.5492E-09 


10.4109 


1.03E-08 


4.404E+04 


0.07 


41.23 


1150 


c 


5.18E+00 


1.69348 


1.4995E-09 


10.4125 


9.75E-09 


4.350E+04 


0.07 


40.87 


1151 


c 


4.89E+00 


1.70032 


1.4498E-09 


10.4142 


9.17E-09 


4.294E+04 


0.07 


40.47 


1152 


c 


4.61E+00 


1.70654 


1.4000E-09 


10.4159 


8.59E-09 


4.236E+04 


0.08 


40.03 


1153 


c 


4.32E+00 


1.71194 


1.3475E-09 


10.4177 


8.02E-09 


4.171E+04 


0.08 


39.52 


1154 


c 


4.00E+00 


1.71555 


1.2921E-09 


10.4198 


7.38E-09 


4.100E-I-04 


0.08 


38.92 


1155 


c 


3.68E+00 


1.71615 


1.2488E-09 


10.4219 


6.74E-09 


4.042E+04 


0.09 


38.41 


1156 


c 


3.51E+00 


1.71368 


1.2020E-09 


10.4231 


6.38E-09 


3.977E+04 


0.09 


37.80 


1157 


c 


3.16E+00 


1.70637 


1.1546E-09 


10.4257 


5.67E-09 


3.909E+04 


0.09 


37.14 


1158 


c 


2.98E+00 


1.69773 


1.1208E-09 


10.4270 


5.31E-09 


3.859E+04 


0.10 


36.66 


1159 


c 


2.79E+00 


1.68575 


1.0858E-09 


10.4285 


4.92E-09 


3.806E-I-04 


0.09 


36.16 


1160 


c 


2.61E+00 


1.67097 


1.0512E-09 


10.4300 


4.55E-09 


3.754E+04 


0.13 


35.67 


1161 


c 


2.43E+00 


1.65364 


1.0168E-09 


10.4315 


4.18E-09 


3.701E+04 


0.12 


35.21 


1162 


c 


2.27E+00 


1.63409 


9.8245E-10 


10.4330 


3.84E-09 


3.649E+04 


0.11 


34.79 


1163 


c 


2.11E+00 


1.61186 


9.4679E-10 


10.4345 


3.49E-09 


3.596E+04 


0.13 


34.41 


1164 


c 


1.95E+00 


1.59444 


9.2036E-10 


10.4361 


3.15E-09 


3.556E-I-04 


0.14 


34.16 


1165 


c 


1.88E+00 


1.58330 


9.0408E-10 


10.4367 


3.00E-09 


3.532E+04 


0.16 


34.02 


1166 


c 


1.81E+00 


1.56475 


8.7896E-10 


10.4374 


2.85E-09 


3.496E+04 


0.16 


33.84 


1167 


c 


1.68E+00 


1.54394 


8.5291E-10 


10.4389 


2.56E-09 


3.458E+04 


0.16 


33.68 


1168 


c 


1.61E-t-00 


1.52883 


8.3498E-10 


10.4396 


2.41E-09 


3.433E-I-04 


0.17 


33.59 


1169 


c 


1.55E+00 


1.51274 


8.1655E-10 


10.4404 


2.26E-09 


3.408E+04 


0.19 


33.52 


1170 


c 


1.48E+00 


1.48977 


7.9113E-10 


10.4412 


2.11E-09 


3.373E+04 


0.19 


33.43 


1171 


R 


1.38E+00 


1.46832 


7.6805E-10 


10.4425 


1.87E-09 


3.342E+04 


0.00 


33.38 


1172 


R 


1.34E+00 


1.44895 


7.4758E-10 


10.4431 


1.77E-09 


3.316E+04 


0.00 


33.34 


1173 


R 


1.25E+00 


1.42850 


7.2624E-10 


10.4442 


1.57E-09 


3.289E+04 


0.00 


33.32 


1174 


R 


1.21E+00 


1.41409 


7.1132E-10 


10.4448 


1.47E-09 


3.270E+04 


0.00 


33.31 


1175 


R 


1.17E+00 


1.39428 


6.9089E-10 


10.4454 


1.37E-09 


3.245E+04 


0.00 


33.31 


1176 


R 


l.lOE+00 


1.36990 


6.6582E-10 


10.4465 


1.20E-09 


3.215E+04 


0.00 


33.32 


1177 


R 


1.05E+00 


1.34597 


6.4121E-10 


10.4474 


1.05E-09 


3.187E+04 


0.00 


33.35 


1178 


R 


9.87E-01 


1.32442 


6.1900E-10 


10.4485 


8.97E-10 


3.162E+04 


0.00 


33.38 


1179 


R 


9.44E-01 


1.30661 


6.0086E-10 


10.4493 


7.84E-10 


3.141E+04 


0.00 


33.42 


1180 


R 


9.07E-01 


1.28575 


5.8007E-10 


10.4500 


6.86E-10 


3.119E+04 


0.00 


33.47 


1181 


R 


8.57E-01 


1.26672 


5.6146E-10 


10.4510 


5.48E-10 


3.099E+04 


0.00 


33.52 


1182 


R 


8.32E-01 


1.25106 


5.4637E-10 


10.4515 


4.81E-10 


3.083E-I-04 


0.00 


33.57 


1183 


R 


7.98E-01 


1.23145 


5.2773E-10 


10.4522 


3.85E-10 


3.064E+04 


0.00 


33.64 


1184 


R 


7.63E-01 


1.20815 


5.0587E-10 


10.4531 


2.84E-10 


3.042E+04 


0.00 


33.74 
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Grid point 


STAT 


T 


X 


P 




M. - M, [g] 


T[K] 


Vc [km s"^] 


Cs [km s ^] 


1185 


R 


7.21E-01 


1.18955 


4.8863E-10 


10.4541 


1.60E-10 


3.025E+04 


0.00 


33.84 


1186 


R 


7.04E-01 


1.17182 


4.7235E-10 


10.4545 


l.llE-10 


3.010E+04 


0.00 


33.93 


1187 


R 


6.67E-01 


1.15923 


4.6089E-10 


10.4556 


O.OOE+00 


2.999E+04 


0.00 


34.01 



